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Abstract 

In stochastic systems with weak noise, the logarithm of the stationary distribution becomes propor¬ 
tional to a large deviation rate fnnction called the quasi-potential. The quasi-potential, and its character¬ 
ization through a variational problem, lies at the core of the Freidlin-Wentzell large deviations theory |32| . 
In many interacting particle systems, the particle density is described by fluctuating hydrodynamics gov¬ 
erned by Macroscopic Fluctuation Theory [^, which formally hts within Freidlin-Wentzell’s framework 
with a weak noise proportional to 1/\/iV, where N is the number of particles. The quasi-potential then 
appears as a natural generalization of the equilibrium free energy to non-equilibrium particle systems. A 
key physical and practical issne is to actually compute quasi-potentials from their variational characteri¬ 
zation for non-equilibrium systems for which detailed balance does not hold. We discuss how to perform 
such a computation perturbatively in an external parameter A, starting from a known quasi-potential for 
A = 0. In a general setup, explicit iterative formulae for all terms of the power-series expansion of the 
quasi-potential are given for the first time. The key point is a proof of solvability conditions that assure 
the existence of the perturbation expansion to all orders. We apply the perturbative approach to diffusive 
particles interacting through a mean-held potential. For such systems, the variational characterization of 
the quasi-potential was proven by Dawson and Gartner [181119| . Our perturbative analysis provides new 
explicit results about the quasi-potential and about huctuations of one-particle observables in a simple 
example of mean held diffusions: the Shinomoto-Kuramoto model of coupled rotators m- This is one of 
few systems for which non-equilibrium free energies can be computed and analyzed in an effective way, 
at least perturbatively. 


1 Introduction 


Large deviations theory studies the exponential decay of probabilities of large fluctuations in stochastic 
systems. Such probabilities are important in many fields, including physics, statistics, finance or engineering, 
as they often yield valuable information about extreme events far from the most probable state or trajectory 
of the system. [HU UHl IMj • 

Weak noise large deviation theory has been developed in the '70s by Freidlin and Wentzell [32] in a 
mathematical framework and by Graham and collaborators [36] with physicists’ perspective. It concerns 
the study of large fluctuations in dynamical systems subject to weak random noise. In this framework, the 
stationary probability to observe some state x of the system obeys a large deviation asymptotics 


-Poo (a;) - exp 



( 1 . 1 ) 
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where e denotes the noise strength squared. The function F is called the quasi-potential. It generalizes the 
notion of (free) energy to general finite-dimensional systems where detailed balance does not hold.but noise is 
weak. If known, the quasi-potential permits to calculate, at leading order in e, important statistical quantities 
such as the probability to observe an arbitrary large fluctuation of the system or the mean residence time 
that the system spends close a metastable state. 

In the last 15 years, Jona-Lasinio and collaborators developed a framework to study large deviations of 
macroscopic quantities (like particle density or current) in a class of many-body systems. Their approach, 
known as the Macroscopic Fluctuation Theory, was mainly applied to stochastic lattice gases, see [5] for a 
recent review. Without stress on mathematical rigor, the MFT can be understood as a generalization of the 
Freidlin-Wentzell theory to the fluctuating diffusive hydrodynamics where a weak noise is added to non-linear 
diffusion equations. Indeed, it is easily obtained by employing a saddle-point approximation in a path integral 
formalism. Even if this approach is only formal because a mathematical meaning of non-linear fluctuating 
hydrodynamics is still lacking, the results obtained are in complete agreement with the ones obtained through 
rigorous probabilistic methods in all cases where a comparison is possible. In this paper, a further example 
of such an agreement will be discussed in chapter 

The Freidlin-Wentzell theory provides a variational characterization of the quasi-potential through dynam¬ 
ical large deviations. However, its explicit computation is typically very difficult and this strongly restricts 
the practical applicability of the theory. This is true both for finite-dimensional and infinite-dimensional 
systems. Only in few cases the quasi-potential can be evaluated analytically, the most important one being, 
of course, the case where detailed balance holds. 

The main focus of the present paper is to answer a very natural but still open question. Let us consider a 
system depending on a control parameter A and suppose that we are able to calculate the quasi-potential for 
A = 0. Can we build a perturbation theory to calculate the quasi-potential for small but finite A? Surprisingly, 
only few works in the literature discuss this question [5^1 US Ell ESI Ell ESI EH EH EH EH • All of them focus 
on specific examples, perform only the l®*-order analysis and most of them consider only finite-dimensional 
dynamical systems (see however m for an example where an infinite dimensional system is analyzed, and the 
case of transition rates between two basins of attraction is studied). We obtain here a precise answer to the 
above question, giving an explicit iterative formula for computing each order in the power series expansion 
of the quasi-potential in A. 

We first analyze the perturbation theory for finite dimensional systems, where the discussion can be made 
quite precise. It is well known that, given the quasi-potential function, one can deduce a simple l®*-order 
equation for the instanton (the corresponding variational problem minimizer). The converse is also obviously 
true: the values of the quasi-potential (the minima) can be easily computed from the minimizers. It is 
however often difficult to compute either the instantons or the quasi-potential without the knowledge of the 
other. We show that this loop can be broken in a perturbative setting: the quasi-potential at any order 
may be computed just from the knowledge of the instanton dynamics of the unperturbed problem and the 
quasi-potential at previous orders. This gives a very natural iterative procedure. We also explain that an 
equivalent simple recursive scheme appears when starting from the perturbative expansion for solutions of 
the Hamilton-Jacobi equation. A key point is to prove that solvability conditions hold at each order, assuring 
the existence of the perturbative expansion to all orders. We show that such conditions are related to the 
behavior close to the attractor, which gives a simple proof that they are satisfied at all orders. 

In the second part of the paper, we consider a particular class of non-equilibrium many-body systems 
described by the Macroscopic Fluctuation Theory [5], namely diffusive particles interacting through mean- 
field potential and driven by an external non-conservative force. Our approach is based on the fact that for 
such diffusions, it is possible to derive a fluctuating hydrodynamics describing the evolution of the empirical 
density for large but finite number of particles N. This evolution equation was first obtained by Dean in 
and it was thereafter called the Dean equation. Although its mathematical status is uncertain, the Dean 
equation allows to treat the mean-field diffusions, at least at a formal level, as a dynamical system in an 
infinite dimensional space perturbed by a weak noise whose strength is proportional to 1/y/N. In the N ^ oo 
limit, similarly to the law of large numbers, the dynamical system becomes a deterministic equation known 
under the name of the McKean-Vlasov or the Vlasov-Fokker-Planck one, as was proven together with the 
propagation of chaos in [5H EH EZj ■ 

We subsequently move our attention to the large deviations around the N ^ oo behavior. The Dean 
equation is a formal random dynamical system to which we apply the Martin-Siggia-Rose formalism [65]. 
Using the saddle point approximation we then end up with an infinite-dimensional generalization of the 
Freidlin-Wentzell theory. This result has been previously obtained rigorously in the mathematical literature 
by Dawson and Gartner EH EH- The formal approach based on the Dean equation makes an explicit 
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connection with the Macroscopic Fluctuation Theory and adds another class of systems to the ones covered 
by the latter theory. 

Explicit results, original to the best of our knowledge, are discussed in the case of the Shinomoto-Kuramoto 
model, a specific stochastic system describing coupled planar rotators m- This model may be also viewed 
as a non-equilibrium version of the dynamical mean-held XY model. We discuss in detail both its TV —?> oo 
behavior and the large deviations around it. For what concern the fV —> oo limit, we are able to fully describe 
analytically the phase diagram of the Shinomoto-Kuramoto model, deriving self-consistent equations for the 
stationary states that are easily solved numerically. 

The perturbation theory developed in the hrst part of the paper is subsequently applied to the Dean 
equation corresponding to the Shinomoto-Kuramoto model, resulting in an explicit calculation of the quasi¬ 
potential close to the known cases: around the free particle dynamics and in the vicinity of stationary states 
of the N = oo theory. This turns out to be a rather simple numerical task. We show that explicit results 
could be obtained to any order but, for the sake of brevity, we only present, the 1®* order calculations of the 
quasi-potential. 

The paper is organized as follows. After Introduction, we begin with a review of Freidlin-Wentzell theory 
results in Sec.[^ Our aim is to present the known results that can be found in |32j with a physicist perspective, 
but nevertheless being precise on many important points that are typically overlooked in the physics literature. 
In particular, we discuss under which hypothesis the quasi-potential can be seen as a solution of the Hamilton- 
Jacobi equation that forms the basis for the perturbative treatment that is developed in Sec.[^ In the 
context of finite-dimensional systems, the discussion of the perturbation theory for quasi-potential is made 
quite precise. We also show that the Taylor expansion of the quasi-potential close to an attractor of the 
deterministic dynamics is a particular case of our perturbative analysis. We conclude the section by observing 
that close to a codimension-one bifurcation of the deterministic dynamics, fluctuations diverge with the mean- 
field critical exponent. Final Sec.is about iV-body systems and is central to this paper. It concentrates on 
the case of overdamped diffusions with mean-field interaction that are introduced in Sec. |4.1| In Sec. |4.2[ we 
discuss the fluctuating hydrodynamics for such diffusions, deriving the Dean equation making the original 
argument [52] more precise. The McKean-Vlasov equation is described as the limit of the Dean equation 
when N —>■ oo. In Sec. 4.3 we introduce the Shinomoto-Kuramoto model m as a particular example of the 


mean-held diffusions and we analyze its McKean-Vlasov limit and the corresponding phase diagram. Sec. 4.4 


derives formally, starting from the Dean equation, an extension of the Freidlin-Wentzell theory to the case of 
mean-held diffusions, and discusses the inhnite-dimensional version of the Hamilton-Jacobi equation for the 
quasi-potential and, briehy, the large deviations for empirical currents. In Sec. |4.5[ we adapt the perturbative 
approach developed for hnite dimensional systems to the diffusions with mean-held interaction, concentrating 
on the power series in the mean-held coupling constant and the Taylor expansion close to stationary solutions 
of the McKean-Vlasov dynamics. A discussion of explicit results for the Shinomoto-Kuramoto model obtained 
by combining analytical and numerical treatments is presented. Sec. [^summarizes the results and discusses 
the perspectives for future work. Two Appendices contain some additional material. 


2 Freidlin-Wentzell theory: a brief summary 

We consider in this section a hnite dimensional random dynamical system dehned by the Ito stochastic 
diherential equatioiQ 


X = K{x) + V^g{x)rit, (2.1) 

where x G g is a d x m matrix and rjt G is a vector of white in time Gaussian noises ry) with zero 
mean and covariance 


= ( 2 . 2 ) 

Throughout the paper it is supposed that K and g are smooth and the Ito stochastic convention is employed, 
if not stated otherwise. We use the notation Q{x) = g{x) g'^{x) with the superscript T indicating the matrix 
transposition and we demand that Q{x) be positive dehnite for all x. Generalizations to the case in which Q 
is semi-positive dehnite are possible but we do not consider this situation here. 

^Here and below, we use physicists’ notation for stochastic equations rather than mathematicians’ one with differentials. 
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We assume that the stochastic process solving Eq. (2.1) has a unique invariant measure. Then its density 
Poo {x) is smooth and solves the stationary Fokker-Planck equation 




i=i 


= 0 . 


(2.3) 


We are interested in the behavior of the above stochastic dynamical system in the small noise limit e <S; 1. In 
this limit, the stationary measure obeys the large deviation principle (1.1) [32j . where the symbol x stands 
for the asymptotic logarithmic equivalence: 


F(x) = — lim e In Pqo (x) . 

e-fO 


(2.4) 


The rate function F{x) is called quasi-potential associated to the random dynamical system (2.1). It is of 


central importance to this paper that P, under suitable hypothesis that will be specified in Sec. |2.3| is the 
unique solution of the Hamilton-Jacobi equation 

VF-[Q-VF + K]{x)=0. (2.5) 

At an informal level, this can be guessed by looking for solutions of the stationary Fokker-Planck equation 


(2.3) obeying Ansatz 


Poo (x) ~ exp 


Fix) 


■Z(x) 


( 2 . 6 ) 


and retaining only terms of order 1/e. The reader may consult |36| where also the lower-order equation for 
Z(x) is derived. 

The present introductory section is devoted to an informal discussion of the Freidlin-Wentzell theory, 
with the particular emphasis on how the quasi-potential can be described from different points of view. This 
knowledge will be used in the subsequent section. We refer the reader to |32] for a mathematical presentation, 
to Chapter 6 of |36j for a treatment more oriented towards the physics community and to [48] for a general 
review oriented to computational aspects. In Sec. |2.1[ we discuss dynamical large deviations, that is, the 
probability that a solution of Eq. (2.1) is close to a given path in the limit e —>■ 0. This permits to define 


quasi-potential Fa relative to an attractor A of the deterministic dynamics x = K(x). In Sec. 2.2 we describe 


the properties of solutions of the Hamilton-Jacobi equation (2.5) that we denote F^j- The connection between 


Fffj, Fa, and F and the conditions under which they coincide are discussed in Sec. |2.3[ concentrating on the 
case of attractive points. A proof of the local existence and uniqueness of F^j around a stable fixed point of 


the deterministic dynamics x = K{r), and of its local regularity, is sketched in Sec. 2.4 


2.1 Freidlin-Wentzell action and quasi-potential relative to an attractor 

Freidlin and Wentzell considered the probability for a trajectory of the stochastic process x{-) defined by 
Eq. (2.1) to be arbitrarily close to a given continuous path x{-) on the time interval They showed 

rigorously that 


sup x(t) — x{t)\ < S 

= lim lim sup e In P 

sup \x{t) — x{t)\ < S 





lim lim inf e In P 
5^0 e^O 


where | • | denotes the norm of a vector in and 

'f 


= At) - Ki^it))] ■ Qixit)) 1 [xit) - Kixit))] dt 


( 2 . 8 ) 


is the so-called Freidlin-Wentzell action functional, see Theorem 2.3 in Chapter 3 of [32]. Observe that such 
a functional vanishes on solutions of the deterministic equation 

X = Kix). (2.9) 

It measures the difficulty for a trajectory to deviate, due to a weak noise, from the deterministic behavior. 
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The above result may be elucidated in the framework of the formal path-integral approach going back to 
Onsager and Machlup [SI [771 ESI na EH- In this approach, the transition probability from a state xi at time 
t = ti to a state Xf at time t = tf is written as the path integral 


\x-i , t j, X ^, t j] 


1 


f 

r A[a:(-)]1 

/ exp 

e 


V[x{-)], 


( 2 . 10 ) 


where the functional integration is restricted to the paths 9 t i— >■ x{t) such that x{ti) = Xi and 


x{tf) = Xf. Za is the normalization factor. The exponential factor exp 
role of the probability density in the space of paths. 


■Afx(-)] 


plays then for small e the 


Let us now consider an attractor A of the deterministic dynamical system (2.9) and a point xq S A. The 
quasi-potential relative to xq is defined as 


= min -^[^(•)] > 

{x(-) I £( —oo)=xo , :e(0)=x} 


( 2 . 11 ) 


where the minimum is over all absolutely continuous paths starting from xq at time t = —oo and ending in 
X at time t = 0. The choice of the time interval [—oo,0] is arbitrary and any semi-infinite interval [—oo,t/] 
would give the same result. It is easy to see that does not depend on the choice oi xq & A and is constant 
on attractors. Indeed, given two points on the attractor, one can always find trajectories starting and ending 
arbitrarily close to them on which the action is arbitrarily small. As an example, the reader may consider a 
limit cycle where any two points of the attractor can be connected by a solution of the deterministic evolution 
X = K{x). For an attractor A, we shall denote by Fa the quasi-potential relative to any xq G A. For the 
sake of simplicity, we shall mainly consider in this paper single point attractors referring the reader to [36j 
and references therein for explicit calculations of quasi-potentials with respect to non-trivial attractors, for 
example, limit cycles or the Lorentz attractor. 

The quasi-potential F can be built from the quasi-potentials relative to the attractors of i = K{x). In 
the case where a fixed point x is the only attractor of the deterministic dynamics x = K(x) for any initial 
condition, F^ coincides with the quasi-potential F defined through the invariant measure in (2.4). Formally, 
this may be understood by observing that P[a;i, ti;xf,tf] is actually independent of Xi in the tf ^ oo limit and 
thus Pooix) =¥[xi,ti', X, oo]. Then, one obtains the stated result by applying the saddle point approximation 
to the right hand side of (2.10). The reader may consult Theorem 4.3 in Chapter 4 of [32] for the rigorous 
result. 

If more attractors {Ai}i<i <7 are present then the quasi-potential F may still be constructed once the 
quasi-potentials relative to each attractor Fh. is known. One has 


F(x) = min [Fa^ (x) + Ci) — min Q 


( 2 . 12 ) 


where constants Ci describe the “height” of each attractor Ai. More precisely, one has to consider the i- 
graphs G{i) on the set {Ai,..., A/} of attractors composed of arrows j ^ k such that j ^ i, from every 
Aj ^ Ai starts exactly one arrow and there are no closed cycles. Constants Ci are then defined as 




(2.13) 


where the minimum is over all z-graphs and the sum over all the arrows of an z-graph and FjWAfe) = Fa^ (x) 
for any x G A^. This rigorous result is discussed in detail in Chapter 6 of [32]- Eq. (2.12) balances the 


contributions from different attractors with the use of the invariant measure of the Markov chain with 
transition probabilities describing the passages between different attractors. This point was discussed in p5] 
from physicist’s perspective. In that reference, a computation of constants Ci for few explicit examples was 
also carried out. In the present paper, we concentrate on the calculation of quasi-potentials F^ relative to 
attractive points x and thus we do not enter into further details on how the heights Ci may be practically 
found. 


2.2 Transverse decomposition, fluctuation and relaxation dynamics 

In the previous section we saw that the quasi-potential F can be obtained by solving the variational problem 


given by Eqs. (2.11), (2.12) and (2.13). In this section, we discuss a different approach that will permit to 


obtain F as a solution of the Hamilton-Jacobi equation (2.5) of Sec. 2.3 
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Consider an open set D CR‘^ and its closure D — DVJdD, where dD denotes the boundary of D assumed 
to be smooth. We suppose in this section that the vector field K admits a transverse decomposition in D in 
the following sense: there exists a smooth function Fhj{x) such that 

K{x) = -{QVFhj){x) + G{x) (2.14) 


and 


VFhj{x) ■ G{x) = 0 


(2.15) 


for all cc S D. The existence of a transverse decomposition is equivalent to demanding that Fhj solves 
the Hamilton-Jacobi equation (2.5). Indeed, if K admits a transverse decomposition, then VFhj ■ K = 


—V Ffij ■ Q V Ffij. Conversely, we can define G = K + Fuj and from Eq. (2.5) we obtain the transversality 


condition. The term QVFhj may be viewed as the gradient of Fhj in the Riemannian metric defined by the 
matrices Q{x)~^ and then Eq. (2.15) states its orthogonality to G with respect to the corresponding scalar 


product of the vector fields. Such an interpretation is often employed in the mathematical literature but we 
shall not pursue it here. 

The deterministic dynamics (2.9) with itr(a:) = —{QWFhj){x) + G{x) is called relaxation dynamics. The 


reason for this name is that as the Freidlin-Wentzell action vanishes on its trajectories, they are in the small 
noise limit the most probable trajectories relaxing to an attractor. The presence of a transverse decomposition 
permits to define the so-called fluctuation or instanton dynamics 


X = 2{QWFhj){x) + K{x) = {QWFhj){x) + G{x) = Kr{x) 


(2.16) 


which plays a fundamental role in what follows. We shall also be interested in the Freidlin-Wentzell action 
functional corresponding to the stochastic dynamics 


X = Kr + \f2eg{x) rjt , 


(2.17) 


which is 


A[a;(-)] = I 


[i(t) — Kr.{x{t))] ■ Q{x{t)) ^ [x{t) — Kr{x{t))] dt. 


(2.18) 


The fluctuation dynamics is connected to time-reversal of the stochastic equation (2.1). To understand 
this point, let us consider the diffusion process defined by (2.1) in the time interval [0,T]. It was shown in 
[47] that its time-reversal corresponds to the stochastic dynamics 


X = —K{x) + V^g{x) r]t , 


(2.19) 



We now discuss some properties of Fhj in connection to relaxation and fluctuation dynamics. The 
stationary points of Fhj correspond to zeros of the vector field K. Indeed 


{K ■ Q-^K)ix) = (VFhj ■ Q VFhjKx) + {G ■ Q-^G){x) > 0, 


( 2 . 20 ) 


as it can be directly proven by inserting Eq. (2.14) into the left hand side and using the transverse decom¬ 


position. Then, if x is such that K{x) = 0, the above expressions imply that VFhj{x) = 0 = G{x). The 
converse is also true if we suppose that the Hessian matrix of Fhj, denoted by (VVFffj), is invertible at x. 
Indeed, by taking the gradient of the transversality condition and evaluating it at x, we infer that 


(VVFjfjG)(x) =0. (2.21) 

By invertibility of the Hessian matrix, this implies that G(x) = 0 and hence K{x) = 0. Moreover, Fhj is a 
Lyapunov function for both the relaxation dynamics and the time-reverse of the fluctuation dynamics 


X = —Kr{x). 


( 2 . 22 ) 
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Streamlines of 

T = K(x) 

Streamlines of 

X = -Kr{x) 


Figure 1: Example showing a simple case where the basins of attraction of the relaxation and of the time- 

2 4 2 

reversal of the fluctuation dynamics do not coincide. Here, Fh,j{xi,X 2 ) = + Q = I, G = 

idFHj/dx 2 , —dFnj/dxi), where I is the 2x2 unit matrix. The dashed blue arrows indicate the direction of 
the relaxation dynamics x = K{x) = —Q\/Fhj{x) + G{x) and the continuous red lines the direction of the 
time-reversal of the fluctuation dynamics x = —Kr{x) = —Q'^Fhj{x) — G(x). 


Indeed, if i: = K(x) then 


dFffj(x) 

dt 


= -i^FHj-QVFHj)ix)<0, 


(2.23) 


and analogously for Eq. (2.22). From the above properties, we conclude that the relaxation and fluctuation 


dynamics have the same stationary points but attractors are transformed into repellers and vice-versa. Thus 
the relaxation and the time-reversal of the fluctuation dynamics have the same attractors and Fjjj as a 
Lyapunov function. 

One could be led to an incorrect conclusion that basins of attraction of those two dynamics are the 
same which, however, is not true because the transverse parts of dynamics can push the systems to different 
attractors. This can be checked, for example, in a very simple two-dimensional bistable system, see Figure[^ 


2.3 Quasi-Potential as a solution of the Hamilton-Jacobi equation 

Let us suppose that Fjjj is a smooth solution of the Hamilton-Jacobi equation 


^Fhj ■ [Q ^Fhj + K]{x)=0, (2.24) 

m D = D U dD C containing in the interior a fixed point x such that Fhj{x) = 0 and Fhj{x) > 0, 
S/Fnjix) ^ 0 for X G D, X X. We shall assume that D is bounded and connected. We want to understand 
the relation between F^ and Fhj, see Theorem 3.1 in Chapter 4 of [35] for more details. 

For any path [—oo, 0] 9 1 1 —>■ x(t) € D subjected to the boundary conditions x{—oo) = x and a;(0) = x, a 
simple manipulation shows that 


A[x{-)] = A[i(-)] + [ ■ ^FHj{x{t)) dt = A[^(')] + Fhj{x) > Fhj{x) (2.25) 

and we used the transversality condition. The inequality 

.A[x(-)] > Fhj^x) (2.26) 

still holds if we allow the trajectory x{-) to go out of the closed set D provided that x G Dhj, the closure of 
the open set 


where Ar[xi-)] was defined in (2.18 


F>hj = {x G D \ Fhj{x) < min Fnjiy)}- 

yGoD 


(2.27) 


Indeed, writing the inequality ( 2.25| ) for the trajectory x(-) restricted to [—oo,r], where r is the first exit 
time from D, we obtain the lower bounds 


-4.[x(-)] > Fhj{x{t)) > Fhj{x) 


(2.28) 
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with the second one resulting from the condition x € Dhj- Action Ar[i(-)] attains its minimum equal to 
zero on the trajectory of the fluctuation dynamics x = Kr{x). It is easy to see that for x € Dhj there 
exists a unique such trajectory [—oo,0] 9 1 1— >■ x(t,x) for which x{—oo,x) = x and i(0,a;) = x. Besides, such 
trajectory lies entirely in Dhj- This follows from the fact, that on each solution of the fluctuation dynamics 
that ends at x, function Fhj decreases backward in time and such a solution may be infinitely extended in 


negative time direction until it reaches x dXt = —oo. The resulting trajectory saturates the inequality (2.25). 
One infers that for x S Dhj-, 


Fhj{x) = 


{x(-) I £( —oo)=ai, ;r(0)=a:} 


A[x{-)] = Fs{x), 


(2.29) 


see (2.11). Dhj is again a bounded open connected neighborhood of x. Its boundary ODhj is composed of 


points oi X G D for which Fhj{x) = min^ga^) Fnjiy)- It is smooth since '^Fhj 0 at such points. 

We have just obtained two results. First, a smooth solution of the Hamilton-Jacobi equation (2.5) on 
with the properties stated at the beginning of the section, coincides with F^ on Dhj where Dhj is the 
sub-domain of D containing x defined by (2.27). Second, on Dhj, 


Fhj{x) = Fs{x) = A[x{- ,x )], 


(2.30) 


where [—oo,0] 9 t i—)■ x(t,x) G Dhj is the solution to the instanton dynamics, i.e. the unique trajectory 
joining a; to x and satisfying Eq. (2.16). 


2.4 Hamiltonian picture 

Above we have assumed the local existence of a smooth solution of the Hamilton-Jacobi equation around an 
attractive point x of the relaxation dynamics (2.9). The linearization of such dynamics around x has the 
form 


X = Ax , 


(2.31) 


where A = (VAr(x))^ is the matrix with entries = Vj Ar®(x). We shall assume that all eigenvalues of A 
have negative real parts, which ensures the exponential convergence to x in the vicinity of the attractor. Such 
attractive points will be called non-degenerate. We shall sketch below a proof of the fact that around non¬ 
degenerate attractive point x there exists a unique local smooth solution of the Hamilton-Jacobi equation. 
The argument we present here is based on the analysis of the dynamics of extremal trajectories of the 
Freidlin-Wentzell action functional, see [2T] for a more global discussion. 

Let us start by considering the Hamiltonian iJ(x,p) related by the Legendre transform to the Lagrangian 


L{x,x) = j{x — K{x)) ■ Q{x) ^{x — K{x)) 


appearing in the Freidlin-Wentzell action (2.8). One has 


F{{x,p) = p • X — L{x,x) = p ■ Q{x)p + p ■ K(x) 


(2.32) 


(2.33) 


for p = ^Q(x) ^(x — K{x)). The Euler-Lagrange equations for the extremal trajectories of the Freidlin- 
Wentzell action correspond to the Hamilton equations 


X = VpiJ(x,p) = 2Q{x)p -I- K{x ), p = -VxH{x,p) = -Vxip ■ Qix)p + p- K{x)). 


(2.34) 


The dynamical system (2.34) in the phase space possesses a hyperbolic fixed point with (x,p) = (x, 0). 

Indeed, the right hand sides of Eqs. (2.34) vanish at this point and the linearization of (2.34) around (x,0) 
has the form 


Phase space 


x = 2Q{x)p + Ax, p=—A^p. (2.35) 

may be split in a unique way into as a direct sum 14 014 of the stable and unstable invariant 


subspaces of the linearized flow (2.35), 


14 = {(x, 0) I X e , 14 = {(a;, Bx) I X G , 


(2.36) 
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Figure 2: Hamiltonian flow around (a;, 0) with the stable and unstable manifolds 


where 


OO 


(2.37) 


is the unique positive definite matrix satisfying the relation 

BA + A^B = -2B Q{x)B . 

On Vs and the linearized flow reduces to 

{x, 0) = {Ax, 0), (i, Bx) = {—B~^A^Bx, —A'^Bx) 


(2.38) 


(2.39) 


and is, respectively, exponentially contracting and exponentially expanding (the eigenvalues of —B~^AB are 
the negatives of those of A and have positive real parts). The subspaces Vs and 14 are Lagrangian, i.e. the 
symplectic form uj = dp ■ dx vanishes when restricted to each of them. It follows from a general theory of 
hyperbolic fixed points, see e.g. [UliO], that in the vicinity of (a;,0) there exist unique stable and unstable 
smooth d-dimensional submanifolds Ms and Mu composed of close points tending to {x, 0) under, respectively, 
flow (2.34) and its time-reversal. Mg is a local piece of 14 around x and Mu has 14 as the tangent space 
at (5,0), see Figure]^ Besides, if K and Q are (real) analytic then the submanifolds Mg and Mu are also 
analytic, see Theorem 7.1 in [49] . 

Both Mg and Mu are Lagrangian submanifolds of the phase space (i.e. 2-form w vanishes when restricted 
to their tangent subspaces). This follows from the Hamiltonian nature of the flow (2.34) which preserves 
UJ. Besides the Hamiltonian H (conserved by the flow (2.34)) has to vanish both on Mg and on since it 
vanishes at (x, 0). We may now define for a: in a small ball D around x 


r 

F{x) = p • dx 
J X 


(2.40) 


with the result independent of the integration path {x{t),p{t)) in such that x{t) lies in D. Then M„ is 
given locally by the equation p{x) = \/F{x), where {x,p{x)) € Mu- Clearly, F is a smooth function on D 
(which is analytic if K and Q are). The vanishing of the Hamiltonian iJ on Mu implies now that F satisfies 


the Hamilton-Jacobi equation (2.5) so that we may set Fhj = F. Note that for such Fhj, the Hamiltonian 


dynamics on the stable manifold Mg projects to the position space to the relaxation dynamics and the one 
on the unstable manifold Mu to the fluctuation dynamics. Fjjj and VFfjj vanish at x and the Hessian 

Fhj{x) = B is positive definite. 

Conversely, if Fhj is a local solution of the Hamilton-Jacobi equation around x with the latter properties 
then, for a small ball D around x, the sets {(a;, 0) | x S D} and {(x, \7Fhj{x)) \ x G D} form, respectively, the 
local stable and unstable manifolds of the fixed point (x, 0) of the Hamiltonian flow (2.34) so that Fhj = F. 

The above argument shows also the local existence of a transverse decomposition (2.14) and (2.15) of the 
vector field K around its non-degenerate stable zeros x. The decomposition has positive definite Hessian 
WFhj{x) and is uniquely determined by this property. 
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The results discussed here, together with the local equality Fhj = T®, show that the quasi-potential 
relative to x is smooth (or analytic for K and Q analytic) in the vicinity of a non-degenerate attractive fixed 
point X of the deterministic dynamics (2.9). It is however well known that does not have to be smooth 
everywhere, see e.g. [HI [371 [S3] and references therein. Such non-smoothness occurs if the unstable 


manifold of the fixed point (a:, 0) of the Hamiltonian dynamics (2.34) has tangent vectors perpendicular to Vg 


leading to the caustics in its projection on 14, see Figurej^ In recent literature, such situations were connected 
to the so called Lagrangian phase transitions (Sj. Non-smooth quasi-potentials may be treated using viscous 
solutions of Hamilton-Jacobi equations [laisi]. We leave such situations to a future investigation. 


3 Perturbative expansions of quasi-potentials 

Let us consider a finite-dimensional stochastic dynamics depending smoothly on a real external parameter A, 

X = K^{x) + g^ix) r]t . (3.1) 

We shall use the notation for the noise covariance and F^ for the quasi-potential relative to 

the attractor x^ of the dynamics x = K^{x) (dropping the subscript indicating the attractor). For simplicity, 
we only consider fixed points as attractors, even if all the results in this section can be generalized to other 
kinds of attractors. 

An explicit calculation of the quasi-potentials F^ for the dynamics of the form (3.1) is usually an impossible 


task. However, it is feasible in particular situations, as for example in the cases when (3.1) respects detailed 


balance. It is thus very natural to ask the following question. Supposing that we are able to calculate the 
quasi-potential for a given value of A, say A = 0, can we perturbatively calculate F^ for small A? Some works 
are available in the literature containing the first-order analysis [39l|40l|4ll|38l|5ll|33l|62] and concentrating 
mainly on specific examples. In |80| . the first-order theory is presented in a general fashion and, recently, a 
rigorous 1®* order analysis has been obtained in m- We extend here the approach of [50] to any order. 

Before discussing in details how the perturbative approach is built in the rest of the section, we first 
summarize the main ideas. The strategy is to consider the quasi-potential as a solution to the Hamilton- 
Jacobi equation 


• [Q^VF^ -f K^] (x) = 0 , 


(3.2) 


which holds under the hypotheses discussed in Sec. 2.3 In particular, we shall assume that the vector field F° 
has a non-degenerate stable zero From the Implicit Function Theorem, it follows that for sufficiently small 
|A| there exists a smooth family x^ of non-degenerate zeros of K^. By a slight modification of the arguments 
in Secs. |2.3] and |2.4[ invoking the dependence on a parameter of the unstable manifold of a hyperbolic fixed 
point, we infer that there exists a smooth family of solutions F^j(x) of the Hamilton-Jacobi equation (3.2) 
defined in a neighborhood of A = 0 and x° and that those solutions coincide with the quasi-potentials 
F'*' for the stochastic dynamics (3.1) relative to the attractors x^ on the set for a sufficiently small 


neighborhood D of x°. D%j is given by Eq. (2.27) for A = 0. Besides F^(x'^) = 0, VF'^(x'^) = 0 and the 


Hessians VVF^(x'^) are positive definite. For K^{x) and Q^{x) analytic in A and x, both x^ and F^j{x) 
will be analytic for sufficiently small |A|. 

For x^ + y G D^j, we shall expand the function A i—)■ F^{x^ + y), well defined and smooth for sufficiently 
small |A|, into the infinite Taylor series around A = 0, writing 


F^(x^ + y) = ^A"F(”)(y) 


(3.3) 


n—0 


Such a Taylor expansion is asymptotic in the smooth case but has a finite radius of convergence in the analytic 
case. An analogous notation will be used for the expansions of and centered at x^. We suppose that 
F(°^ is explicitly known and attempt to calcul ate the perturbative coefficients F^'^^y). Inserting the above 
expansion into the Hamilton-Jacobi equation (3.2), one obtains a hierarchy of equations expressing F*-" ^ in 
terms of functions F^^'> with k < n so that they may be solved iteratively. We stress that in Eq. (3.3) we 
have moved attractors x^ to the origin. This may seem just a detail but it is important in order to g et a 


y&Kj- 


Recall from Sec. 


2.3 


that 


simple proof that the equations for F^") admit a unique solution for x° - 

x° may be connected to each x G by a unique trajectory of the A = 0 fluctuation dynamics which, after 
the shift by —x°, takes the form 


y = Kl°\y) = 2 (Q^^\y)VF^°\y)) + K^^\y). 


(3.4) 


10 













The solution for e ach F^^'> can be obtained with the method of characteristics by integrating along the 
trajectories of (3.4), see Eq. (3.16) below for the explicit expression for . Such obse rvation is of practical 
importance: to compute the only difficulty is to compute solutions of Eq. (3.4). This might not be 

doable analytically, but it is a simple problem for a numerical treatment. 

From the practical point of v iew it may be more convenient to implement a modified version of perturbative 
expansion where we replace Eq. (3.3) with a Taylor expansion of the function A i— F^{x) at A = 0 for x G D'jjj : 


F^{x) = y A"F(”)(a:). 


(3.5) 


n—0 


We shall see an example where expansion (3.5), that we shall call direct, is simpler to implement than 


Eq. (3.3). The difference between the two expansions comes from the fact that here we do not move the 
A-attractors to the origin. In this case, however, the simplest route to prove that functions may be 


again found iteratively with the meth od of characteristics, is to reduce expansion (3.5) to (3.4), see Sec.|3^ 
Finally we shall show in Sec. 3.3 that, for a suitable choice of and the expansion of Eq. (|3.3|) 


reduces to the Taylor expansion of the quasi-potential F associated to Eq. (2.1) around the attractor x. 


3.1 Expansion centered on attractors of the perturbed dynamics 


This section discusses the perturbative solution of the Hamilton-Jacobi equation ( |3.2[ ) in the form (3.3). We 
iteratively construct functions F*^") obtaining, in such a way, two results: first, we prove that the equations 
obeyed by F^”^ admit a unique solution and, second, we give an explicit recursive formula that may be used 
to practically calculate F^"). 

As already mentioned, we assume that xP is a non-degenerate fixed point of the relaxation dynamics 
at A = 0. This means that K^{xP) = 0 and the eigenvalues of VK^(x^) have negative real parts. The 
assumption ensures that quasi-potential F° relative to x^ (recall that we have dropped the subscript indicating 
the attractor) is given by the solution to the Hamilton-Jacobi equation in the neighborhood F>^j of x^, that 
F°(x°) and VF ^(x^ ) vanish, and that the Hessian VVF°(a;°) is positive definite and F°(x) >0 for a: G D^j, 


X 7 ^ a;°, see Sec. 2.3 


Let us expand the drift and the noise covariance in powers of A after shifting the attractor x^ of the 
perturbed deterministic dynamics to the origin: 


F^(a:^ + y) = y A"F(")(y), + y) = ^ A" Q'") (y). (3.6) 

n—0 n—0 

Several properties that will be used below follow simply. Because x^ is a fixed point for the deterministic 
dynamics, F^"'^(0) = 0. Similarly, since F^{x^) = 0 and VF^(a;'^) = 0, the relations F*-"^(0) = 0 and 
VF*^")(0) = 0 must hold for every n. Moreover, the eigenvalues of VF*^°^(0) have negative real parts and 
VVF(°)(0) and are positive definite. 

Inserting expansions (3.3) and (3.6) into the Hamilton-Jacobi equation ( |3.2[ ), we obtain the power-series 
identity 




n—k 


y vfp 


n—0 k—0 _ /—0 

Upon equating to zero order by order, this gives a hierarchy of relations 


VF(°) • qP'>\/fP^ + kP'> =0 


for 

for 


(y) = 0. 


n = 0, 
n > 0, 


where 


r(P = 2QP^VFP^ 

and is a functional of FP\ ..., fP~P given by 


RiO) 


(3.7) 

(3.8) 

(3.9) 

(3.10) 


n—1 r 


/c=l 


n—k 


5'(")[f(0),...^f(”"^^] = -y \/fP~P ■ (qP'>WFP^ + Rp'>'^ 


1=0 
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(3.11) 


We have arranged Eqs. (3.91 in such a way that appears only on the left hand side. Eq. (3.81 is nothing 


else but the Hamilton-Jacobi equation corresponding to the dynamics (3.11 with A = 0, once we have moved 
the attractor to t he o rigin. We assumed that its solution (y) is known. In the following, we prove that 
solutions to Eqs. (3.9) for n > 0 exist and are unique. An explicit formula will permit to obtain given 
F*^^^ for k < n. We start by looking at the properties of the 0*^ order fluctuation dynamics y = (y) when 


the norm \y\ is small. These results will be useful to prove the existence of solutions to Eqs. (3.9). 


3.1.1 Fluctuation dynamics: exponential escape from attractor 


Let us consider Eq. (3.4) describing the 0*^ order fluctuation dynamics after the shift of attractor to the 
origin. 

ye 


From the results of Sec. |2.3| about the trajectories of the fluctuation dynamics it follows that for 
4. — = Do there exists a unique solution [—oo, 0] 9 1 1 -)- y{t, y) of Eq. ( |3.4[ ) such that y{—oo, y) = 0 

and y{0,y) = y. Besides, y{t,y) belongs to Dq for all t. 

We now show that y(t,y) escapes from the attractor y = 0 exponentially fast. Indeed, we can write for 
\y\ small 


gt(2Q<«)B+A) y 


y{t, y) = y + o ^ 

In the above expressions, A and B are defined through the small \y\ expansion of and F^°^: 

= Ay + o{\y \), F^°'>{y) = \y ■ By + o {\y\^) , 


(3.12) 


(3.13) 


i.e. A = (VF(°)(0))^ and B = VVF*^°)(0). One should recall from the properties listed above for and 
F(°) that the eigenvalues of A have negative real parts and B = is positive definite. We have encountered 
matrices A and B already before when studying the Hamiltonian dynamics in Sec. |2.4| They are related by 
the identity 


BA + A^B + 2B Q^°'>B = 0 


(3.14) 


imposed by the 2'^'^ order in y contribution to the 0*^ order Hamilton-Jacobi equation (3.8). Note that 


identity (3.14), which implies that 


B F A = -B-^ A^B , 


(3.15) 


coincides with Eq. (2.38) solved by (2.37) if we replace in the latter Q{x) by Q*^°^(0). Behavior (3.12) is 
dictated by the spectrum of {2Q'^^'^B F A). Indeed, as already noticed before, the eigenvalues of —BA^B have 
positive real parts. Thus the 0*^-order fluctuation dynamics escapes exponentially fast from the attractor. 


This agrees with the analysis of Sec. |2.4| where we showed that the fluctuation dynamics on the position space 
corresponds to the Hamiltonian dynamics on the unstable manifold of the hyperbolic fixed point. 


3.1.2 Iterative solution 


We shall prove that the unique solution of Eqs. (3.9) on Do = — x° that satisfies F(")(0) = 0 is 


F(")(y)= r F('^)[FW,...,F("-i)](y(t,y))dt, 

^ —OO 


(3.16) 


where y{t,y) is the trajectory of the 0**'-order fluctuation dynamics (3.4) joining the origin to y that was 


discussed above. Moreover, the expression (3.16) is well-defined, as we show that the integral in this expression 


is convergent whenever y belongs to Do- Eq. (3.16) give s an iterative solution of (3.9). 

Let us start by proving that if a solution of Eq. (3.9) such that F("’)(0) = 0 exists, it has to have the form 


(3.16). This is easily seen by taking the total time derivative of F^^\y(t^y)) 


|F(")(y(t,y)) = VF^-\y{t,y)) 


(vF(").F(°))(y(t,y)) 
= F(")[F(o),...,F(-i)](y(t,y)), 
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(3.17) 































and by integrating over time. 

In the next step, we shall show that, assuming that the integ ral o n the right hand side of (3.16) converges, 
the latter equation gives a function that solves Eq. (3.9). Indeed, on the fluctuating dynamics 

trajectory. 


F^^Hy{s,y)) = f S^^\---]{y{t,y{s,v)))dt= f ](y(t + s,y)) 

J —oo J —oo 

= r S<^-\---]{y{t,y))dt, 


dt 


(3.18) 


where we used the relation y{t,y{s,y)) = y{t + s,y) that holds because of the uniqueness of the solutions 


of Eq. (3.4). Now, deriving the previous expression with respect to s and evaluating at s = 0, we obtain 


Eq. (3.9). This means that expression (3.16) solves Eq. (3.9) provided that it is well defined. 


To complete the proof, we need to show that the integral appearing in Eq. (3.16) is convergent and vanishes 


at 2 / = 0. To show this, we have to analyze the behavior of the integrand when t —> —oo. This corresponds 


to studying the behavior of the integrand for small \y\. One completes the prove combining Eq. (3.12) with 
the observation that the expression defined in Eq. (3.11) is at least quadratic in y for \y\ small if 
for k < n have the same property. Under this assumption, that is true for F^^\ the integrand on the right 


hand side of Eq. (3.16) converges exponentially to zero when t —>■ —oo so that the time integral converges. 
Besides, it determines function that is at least quadratic in y for small \y\ so that it may be shown 

inductively that starts at worst quadratically. 

Let us summariz e ou r results: we desc ribed an iterative scheme to calculate perturbatively the quasi¬ 
potential F^ of Eq. (3.1) as a power s eries (3.3). Once we know the terms of th at ex pansion for A: < n 
then E^”) may be obtained using Eq. (3.16). In that formula, is defined by Eq. (3.11) and y{t,y) is the 
solution to the 0*^-order fluctuation dynamics (3.4) that starts at the origin at t = —oo and arrives at y 
at time zero. Finally, we proved that this procedure is well defined (E^"^ are finite quantities and depend 
smoothly, or analyticaly in the analytic case, on y € Eg). More generally, one could define this way E^")(y) 
for y belonging to the basin of attraction of the origin for the time-reversal of the fluctuation dynamics ( |3.4[ ) . 

Let us conclude this section observing that the procedure described here may be very easily implemented 
numerically. Indeed, to calcul ate F ^'^'> at all the orders, one only needs to compute the solution to the 
O^ii-order fluctuation dynamics (3.4). The scheme just described gives a powerful practical tool to compute 
quasi-potentials perturbatively. 


3.2 Direct Expansion 


In Sec. |3.1| we have presented an expansion in power of A centered on the attractors of the perturbed dynamics. 
We discuss in this section a direct expansion, which does not depend on the knowledge of the attractors of 
the perturbed dynamics. The main point we want to stress in that case is the appearance of non-trivial 
solvability conditions. 

There are several reason why, in some cases, this new expansion may be simpler or more relevant than the 
one centered on the attractors. The main one is that the attractors of the perturbed dynamics may not be 
known explicitly and should then be co mpu ted themselves by a perturbative expansion. In that case, as the 
definition of and given in Eq. ( |3.6[ ) involves the attractor of the perturbed dynamics x = K^{x), 
then matrices may be non-zero at all orders even if the covariance is independent of A. An example 
of this kind will be encountered in Sec. 14.51 

We assume power series expansions for and 


K^{x) = y \^k^^\x ), Q^{x) = y A”Q(”)(a:). 


n=0 


n=0 


(3.19) 


Observe that = 0 and VE(°^(a;°) = 0. By assumption all the eigenvalues of ViL^°^(x°) have 

negative real part, and is a positive definite matrix. This implies that VVE(°^(a;°) is a positive 

definite matrix. As VK^^^(x^) is invertible, the identity K^(x^) — 0 permits to solve iteratively for the 
coefficients of the Taylor expansion 


= y . 

n=0 


(3.20) 
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for the perturbed attractor, where = x^. 
Inserti ng E g. 


and (3.191 into the Hamilton-Jacobi equation (3.2|, we obtain a hierarchy identical to 


(3.8) and (3.9). The only differences are, of course, that y has to be replaced with x and and 

by and Explicitly, we have 


where 


VF^'^\x)- k^^^Xx) = 

= 2q(°)ve(°) + 


for 

for 


n = 0, 
n 0, 


(3.21) 

(3.22) 

(3.23) 


and 5*-"^ is the functional of F^^\ ..., given by Eq. (3.11) with and replaced by the 

hatted quantities. 

We first remark that, as k^\x'^) = 0, Eq. (3.22) implies 


,§(") [F^^\ ... 


= 0 , 


(3.24) 


which appears as a solvability condition for Eq. ( |3.22 1. It is possible to prove directly this solvability condition 
by induction, however this involves subtle cancellations that are tedious to prove to all orders. In order to 


bypass this proof, we rather argue that we know a priori that the series expansion exists. Then Eq. (3.22) is 


a consequence of the existence of the series expansion and this implies that the solvability condition (3.24) 


is satisfied. The existence of the series expansion follows from the existence of the series expansion around 
the attractors of the perturbed dynamics, discussed in the previous section. Indeed, F^'^\x) may be directly 
found using the expansion from the previous section by comparing order by order both sides of the identity 


OO 

E 

n—Q 


A"F(")(a;) = F^{x) = V \^F^^\x - x^) , 


n=0 


(3.25) 


where on the right hand side one inse rts th e Tayl or exp ansion (3.20). 

We now explain how to solve Eq. (3.22) using ( 3.24| ). Let us consider the 0*^-order fluctuation dynamics 

x = k'f\x), (3.26) 

and its trajectory x{t,x) lying in D'^j such that x{—oo,x) = x^ and i(0,x) = x. With the same argument 


as in Sec. |3.1.2] we can show that 

F(")(a;) = + / F(”)[F(°\...,F("-i)](f(t,a;))dt. 

J —OO 


(3.27) 


has to hold for a solu tion of Eq. (3.22), where are (for the moment arbitrary) constants. This is 
analogous to Eq. (3.16). To obtai n F^'^ '>{x) from this equation, we have to prove the convergence of the 
integral on the right hand side of (3.27) and to fix the constants The first task requires the c ontrol 

of the behavior of ..., Fl”“^->](x) around x = xP. It may be achieved by induction using (3.24) 


and the exponential relaxations of x{t,x) to x^ when t goes to —oo. The second task is easier and may be 
accomplished iteratively since the normalization F^(x^) = 0 leads upon Taylor expending to the relations 
that allow to express = F^'^\x^) b y the values at x^ of functions F^^'> and their derivatives and by the 
coefficients of the Taylor expansion (3.20), all for k < n. Note, however, that the choice of ioi k < n 
in not relevant for the calculation of 
gradients of F^^'i for fc < n that enter 5^ 


except when it comes to the choice of Indeed, these are the 


3.3 Taylor expansion of the quasi-potential around the attractor of the unper¬ 
turbed dynamics 


Consider now the stochastic evolution (2.1) without an external parameter. Let a; be a non-degenerate 


attractive point of the deterministic dynamics (2.9). We shall be interested here in the Taylor expansion of 
the quasi-potential in a neighborhood of x : 


OO 

F{x + y) = Y^ 


(3.28) 


n=2 
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where and are rank-n tensors with the components 


\ / ii ■■■i 


1 

n! 9a:®i- • • 




(3.29) 


and on the right hand side of (3.28) the contraction of all indices is implied. We shall show here that expansion 
(3.28) can be viewed as a particular case of the perturbative expansion considered in Sec. 3.1 Then, all the 
results obtained there may be applied to (3.28) providing a method to calculate tensors 
With the same notation as in (3.28), we introduce the expansions of K and Q around x : 


OO OO 

K{x + y) = Y^ , Qix + y) = Y, . 


(3.30) 


n—1 


n—0 


Since x is a non-degenerate attractive fixed point of the relaxation dynamics, — VK{x) has eigenvalues 

with negative real parts and = ^WF(x) is a positive definite matrix. 

Let us introduce a new A-dependent system defined by 


K^{x + y)=\K{x + Ay) , Q^{x + y) = Q{x + Ay) 


(3.31) 


that reduces for A = 1 to the previous one and depends smoothly on real A. Note that point x is a non¬ 
degenerate stable attractive zero of for all A and that 


F {x + y) = ■^F{x + \y) 


(3.32) 


satisfies the Hamilton-Jacobi equation (3.2) for all A including A = 0. The scaling with A was introduced in 
such a way that the A = 0 case gives a non-trivial contribution. It allows to align the notations to those of 
Sec.|3.1| With our choice, we indeed have 


K^{x -t y) = (y • W)K{x) = Ay , Q°{x + y) = Q{x ), 
F°(x -h y) = i(y • V)^F(x) = \y ■ By . 


(3.33) 

(3.34) 


Moreover 

if(")(y) = , Q(")(y) = (v(”)Q)y( 


F(")(y) = (^v("+2)F)y(’"+2) (3.35) 


in the notation of Eq. (3.6) and (3.3). It is worth stressing that, in this context, the role of the unperturbed 
A = 0 stochastic dynamics is played by the linear dynamic 

y = Ay + g{x) yt (3.36) 


with matrix A given by (3.33), for which the unperturbed quasi-potential is the quadratic approximation of 
F around the attractor. This is not surprising since the stochastic equation (3.36) defines a Gaussian process 
of the Orstein-Uhlenbeck type whose invariant measure is Gaussian with the covariance equal to eB~^, where 
B is the Hessian matrix of F at x, see (3.34). The Hamilton-Jacobi equation reduces for A = 0 to the identity 
(2.38) with the solution given by (2.37). The relaxation and the fluctuation dynamics are (after the shift of 
the attractor to the origin), respectively. 


y = Ay 


and 


y = -B-^A^By. 


(3.37) 


In particular, for A = 0 the trajectories of the fluctuation dynamics that start at t = —oo from the origin 
have a particularly simple form: 


y(t,y)=e-‘^'^^% 


and Eq. (3.16) reduces to the iterative solution 

rO 


F(")(y)= f 

J —OO 




y)dt 


(3.38) 


(3.39) 


for the Taylor coefficients of F at x with given Eq. (3.II). 
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3.4 Codimension-one bifurcations: the critical exponent 


Let us return to a family of dynamical systems parameterized by A as in Eq. (3.1) and let us suppose that 
the deterministic dynamics i = K^{x) has an attractive fixed point for A < Ac that is non-degenerate 
for A < Ac and undergoes a codimension-one bifurcation at A = Ac- Large deviations for normal forms 
corresponding to codimension-one and codimension-two bifurcations were discussed in a series of papers in 
the ’70s and ’80s. We refer to [511 [13 [5S1HI] and to the review jSD] for a detailed analysis. Here we want 
only to remark that the critical exponent in this framework is equal to that of the mean-field theory. 

The assumed scenario implies that \7K^{x^) has a simple real eigenvalue that is negative for A < Ac 
such that 


a 


= 0 


and 


da^ 

~dX 


\=\a 


7^0 


(3.40) 


(the saddle-node bifurcation) or a complex eigenvalue with negative real part for A < Ac, and its conjugate, 
such that 


Rea^‘= = 0, Ima^‘= ^ 0 


and 


da^ 

~dX 


\=\a 


7^0 


(3.41) 


(the Hopf bifurcation). Moreover all the other eigenvalues of VK^{x^) have strictly negative real parts for 

A < Ac. 

Let us consider the covariance 


{c^y^ = 


P^i.x)dx - [J PPL{x)dx^(^ J x^ P^{x)dx^ 


(3.42) 


of the invaria nt m easure of the stochastic dynamics (3.1). If the minimum of the quasi-potential desc ribing 
the behavior (1.1 1 of in the limit of small noise is attained at then the saddle-point analysis of (3.421 
implies that 


poo 

lim -C^ = (P^)-i = 2 / 


(3.43) 


whe re is the Hessian matrix of at and the last equality with = {WK^{x^))'^ follows from 
Eq. (2.37). Hence may be viewed as the small-noise limit of the stationary equal-time connected 

2-point function of the stochastic process solving Eq. (3.1). Let us examine the behavior of the latter limiting 
quantity when A Ac. Any vector y G 


may be decomposed as 
y = G- y^ 

where /3^ G C, is the eigenvector of (A'*')^ with the eigenvalue and y^ G 
subspace of {A^)'^ corresponding to the other eigenvalues so that 


(3.44) 

belongs to the invariant 




(3.45) 


in the vicinity of Ac for some 5 > 0. Then the 2-point function 

y{B^)-^y = -2\P^f{Rea^)-^^-Q^{x^)v^-2Re[{P^f{a^)-'^v^-Q^{x^)v^] -f 0(1) (3.46) 

where 0(1) term is regular when A ^ Ac. Assume that f3^'= y 0 which holds for generic y. Then, when 
A Ac, the first term on the right hand side diverges like (Ac — A)“^ and the second term is regular if the 
purely imaginary eigenvalue a\^ y 0 , and both terms have the same divergence when a\^ = 0 . 

This completes the simple proof that, in the framework of a dynamical systems in the low noise limit 
undergoing a codimension-one bifurcation, critical exponent for two-points correlator is equal to 1 , the mean 
field theory value of the susceptibility exponent 7 . 


4 Mean-field systems and the Shinomoto-Kuramoto model 

We consider in this section systems composed of many diffusive particles, interacting through a mean-field 
type two-body potential and driven out of equilibrium by external forces. Most of our analysis is valid for any 
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model in this class (the two-body interaction could even be non-potential). However, in order to go beyond 
formal results, we discuss in detail the example of Shinomoto-Kuramoto system, a simple 1-d model first 
introduced in [74] , which has recently attracted some attention in the mathematical literature |60l [Ml 0135| 
as well as in the physical one izniiziiiHa. 

It is known since the 60s that the evolution of the empirical measure of N diffusions with mean-field 
interaction is described, for iV —>■ oo, by a non-linear Fokker-Planck equation known in the mathematics 
literature as the McKean-Vlasov equation [55], see also [7niis3 for more modern presentations. We refer 
to this limit as the mean-field behavior of the N particle system. For the Shinomoto-Kuramoto model, 
the mean-field behavior exhibits a rich phase diagram with stationary and periodic phases separated by 
bifurcation lines, as first observed in [75] . 

In [151 [TU] , Dawson and Gartner studied the large deviations of the empirical measure for the mean-field 
diffusions through a generalization of the Freidlin-Wentzell theory to such questions. Extensions to cases 
where quenched disorder is present, as for example in the Kuramoto model [Mim or in mean-field spin 
glasses [33 131 [5], were considered too. Also some large deviation results on a mean-field model for active 
matter were obtained in the physics literature [3]- 

In this chapter, we first discuss how the mathematical results [181 US] can be obtained formally (i.e., 
without mathematical rigor) by writing an effective evolution for the empirical measure in the form of a 
stochastic PDE, called the Dean equation [52|. The noise term in this equation is proportional to I/VN so 
that it vanishes as N —>■ oo. In the latter limit, one recovers the deterministic McKean-Vlasov equation. 


The evolution of the empirical measure is thus formally given by a stochastic partial differential equation 
with weak noise which, at variance with the cases considered in the previous sections, is infinite-dimensional. 
By applying functional integral techniques from field theory (the Martin-Siggia-Rose formalism, see [5S|), we 
can write the infinite-dimensional analogue of the Ereidlin-Wentzell theory. In the physics literature, this 
extension goes under the name of Macroscopic Eluctuation Theory and it has attracted much attention in 
the statistical mechanics community over last years, see 0 for a review. In the present case, we re-obtain 
formally the rigorous results by Dawson and Gartner [181 US] ■ We also discuss a more general large deviation 
result covering current fluctuations. 

We then apply the infinite-dimensional analogue of the perturbative scheme discussed in Ghapter 0 to 
obtain explicit results about the quasi-potential for the mean-field diffusions. In particular, we calculate 
the quasi-potential perturbatively close to the free particle dynamics with no interactions between parti¬ 
cles present. This permits also to compute perturbatively the rate function for the fluctuations of some 
macroscopic observables. Eor the Shinomoto-Kuramoto model, explicit results obtained by implementing a 
numerical algorithm to compute the quasi-potential to the I®* order in the coupling are presented, together 
with an analysis of the fluctuations of magnetization. 


We also discuss how explicit results may be obtained for the Taylor expansion of the quasi-potential around 
a stationary solution of the McKean-Vlasov equation. Within this analysis, it is clear that the variance of 
density fluctuations diverges close to bifurcations when external parameters are changed, a result rigorously 
obtained in [ZQlE]. A numerical algorithm to evaluate explicitly the Taylor expansion is also discussed but 
we have not implemented it. 


The structure of this part of the article is as follows. In Sec. |4.I[ we introduce the class of systems that will 
be considered. Sec.|4.2|derives the Dean equation, formally describing the evolution of the empirical measure 


for large but finite N. Then, in Sect. 4.3 the particular case of the Shinomoto-Kuramoto model is discussed. 
We describe the long-time behavior of solutions of the associated McKean-Vlasov equation. Despite its 
simplicity, the model, which is a kinetic and non-equilibrium version of a mean-field ferromagnet, exhibits 
a rather complex mean-field behavior that we describe concentrating on the results that can be obtained 
analytically or semi-analytically. In Sec. |4.4[ we discuss how a generalization of the Freidlin-Wentzell theory 
may be formally obtained by applying the Martin-Siggia-Rose formalism to the Dean equation. Finally, in 
Sec. |4.5[ the perturbative calculation of the quasi-potential is performed and some explicit results for the 
Shinomoto-Kuramoto model are described. 


4.1 Mean-field diffusions 

Let us consider a system composed of N particles undergoing an over-damped diffusion in or in a torus 
and coupled through a mean-field 2-body potential V{x) = V{—x). The equations of motions defining 
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the stochastic evolution are 


±n = b{Xn) - —'^(yV){Xn-Xm) + \/‘^ksT r]n , (4.1) 

m—1 

where rjnit) are independent white noises with zero average and covariance E 77 )^(t) 77 ^(s) = S^^6n,mS(t — s), 
T > 0 is the temperature, and ks is the Boltzmann constant. 

The quantities which are of central interest for us are the empirical density and the empirical current, 
defined as 


1 X ^ 

PN{t,x) =—'^6{x-Xnit)) (4.2) 

n—1 
1 ^ 

jNit, x) = - Xn{t)) o Xn{t) (4.3) 

n—1 

where ” o ” stands for the product in the Stratonovich convention. It is straightforward to show that the 
following continuity equation holds 


dtPN + V • JAT = 0 


(4.4) 


(this uses the chain rule which imposes the Stratonovich convention in the definition of the empirical cur¬ 
rent). Moreover, substituting the equation of motion (4.1) into the definition of jjv and returning to the Ito 
convention, we obtain 


jN{t,x) = pN{t,x)[b{x) - JV{V * PN){t,x)] - ksTVpN{t,x) + ^ 6(x - a;„(t)) ? 7 „(t), (4.5) 


n—1 


where (V * pn) is the convolution between V and pN{t, •) and the term with VpN was produced by the change 


of the stochastic convention. Observe that Eq. (4.4) with jjv given by (4.5) is not a closed equation for p^ 


because the noise term explicitly depends on the particle positions. What we would like to do, instead, is to 
write a closed evolution equation for into which the particle positions enter only through ppf. 


4.2 Evolution of the empirical density: the Dean equation 


A closed equation for the evolution of the empirical density was obtained from Eqs. (4.4) and (4.5) by Dean 
[22] . Dean’s argument (somewhat brief in the original paper) can be reformulated in the following way. The 


last term on the right hand side of (4.5) may be viewed as a white noise in time with values in vector fields 
on 


- Xn)r]n{t) 


(4.6) 


parameterized by the particle positions {xn)^^i. White noise (4.6) has mean zero and covariance 


‘^^bT ^ ^ d{t - t') 5{x - y) p{x ), 


iV2 


N 


(4.7) 


where the last expression followed by using the distributional identity 5{x — Xn) S{y — Xn) = S{x — y) 5{x — Xn) 
and introducing the particle density p{x) = ~ ^n-)- Consider now another noise. 


2kBT 

N 


y/p(x)Ci.t,x) , 


where ^{t,x) stands for the vector-valued white noise in space and time satisfying 


E C{t,x) = 0, E C{t,x)^^{s,y) = J*-’ d{t- s)6{x - y). 


(4.8) 


(4.9) 
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Random process (4.81 may again be viewed as a white noise in time with values in vector fields on but 
now parameterized by densities p{x). It has a zero mean and covariance 


2kBT 

N 


5'’^ 6{t — t') d{x — y) p{x) 


(4.10) 


that coincides with the one of noise (4.6) if p is related to particle positions as above. Dean proceeded 
identifying the two white noises by writing 


jN{x,t) = jp^{t,x) + y VPN{t, x) C(t, x), 


(4.11) 


upon which the continuity equation (4.4) became a closed stochastic PDE in the space of densities for 
p(t,x) = PN(t,x), 


dtp{t,x) + ■jp{t,x) + ^ 
where jp is the nonlinear functional of p given by 


2kBT 

~7r~ 


p{t,x)^{t,x)^ = 0 , 


jp{t,x) = p{t,x)(b{x) - JJ{\7V){x - y) p{t,y) - ksT^Ip{t,x). 


(4.12) 


(4.13) 


We shall call (4.12) the Dean equation for the empirical density. 

Dean’s substitution is formal in the infinite-dimensional situation involving the space of densities but 
would be legitimate in a finite-dimensional setup. To explain what we mean, let us consider the backward 
Kolmogorov equation describing the evolution of averages for functionals of the empirical density iS[/3jv(t, •)]■ 
Denoting by E the average with respect to the noises pn and applying the Ito calculus, we infer thalj^ 

S[pNit,-)] = J dx -^^[pN{t,-)]'^ ■ (^PNit,x)[-b{x) + JV(y * pN)it,x)] +kBTVpN{t,x)^ 


kBT 


6‘^S 


dxdyj^^^-^^^[pN{t,-)]V,-Vy(5{x-Xn{t))5{y-Xn{t)) 

J ^ (^PN{t,x)[- b{x) + JViy * pN){t,x)] + kBTVpN{t,x) 

[pN{t,-)]Va: ■ Vy(^PNix,t)Six - y)y (4.14) 


kBT 

N 


dx / dy 


5p{x)5p{y) 


where the last expression was obtained proceeding as in (4.7). We thus obtain for E [5[pjv(^, •)]] the evolution 
equation 


^E 5[pjv(t,-)] = E/:5[pAr(t, •)], 


(4.15) 


where C is the generator given by 

^S[p] = jdx^^^V-(^pM{t,x)[-b{x) + JV{y*pN){t,x)\+kBTVpN{t,x)^ 


kBT 

N 


II *■> ^ *> - ">) ^ 


(4.16) 


Observe now that the above evolution for E [5 [pat]] is given by the same backward Kolmogorov equation that 
the one obtained assuming that pN solves the Dean stochastic PDE (4.12) anticipated at the beginning of 


the subsection. In finite dimension, it is not enough to know that Ef{X{t)) = ELf{X{t) for a generator 
of a diffusion process L and each / to deduce that X{t) has the law of the diffusion process satisfying 


^Since the densities are normalized, functional derivatives 5S/5p{x) are defined only up to a constant, but such ambiguities 
drop out in all expressions below where the functional derivatives are integrated against functions with vanishing integral. 
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the corresponding SDE, but this would follow if one showed that f{X{t)) — f{X{0)) — J* Lf{X{s)) ds are 
martingales m- A slight extension of the previous calculation shows that the process 


t 

S[pNit, •)] - 5[pAr(0, •)]] - J CS[pNis, •)] ds 


(4.17) 


possess this property. Nevertheless, since we are in infinite dimension, the derivation of the Dean equation 
remains formal. 

To complicate things further, the mathematical status of the Dean equation is unclear: it is difficult to give 
sense to different terms of the equation, including the noisy one, in a function space that would contain the 


(distributional) empirical densities of coupled diffusions (4.1), not even speaking about a theory of solutions of 
such a stochastic PDE. Although there is a considerable mathematical literature about non-linear stochastic 
PDE’s with the noise that is delta-correlated in time and space, see isomziiis], it does not cover the case of 
the Dean equation. One may hope, however, that it is possible to give a meaning to this equation at least 
in the large deviation regime for large N. Although in what follows we completely avoid such rigorous issues 
and proceed formally, the fact that our results agree with those rigorously obtained by Dawson and Gartner 
[BIS] gives support to this assumption. 

Before proceeding further, let us observe that the Dean equation can be recast in a form similar to that 
of the finite-dimensional SDE considered in the first part of the paper by rewriting it as 


dtp = JC[p\ + \ 1^ V[p], 


where 


K:[p\{t,x) = -V ■ jp{t,x) 

is the drift and 

■n[p]{t,x) = -^/ksT V • (^^/p{x)^{t,x)^ 

is the white noise in time, parameterized by density p, with zero mean and covariance 


(4.18) 

(4.19) 

(4.20) 


E T][p]{t,x)r][p]{s,y) = S{t - s) Q[p]{x,y) 


(4.21) 


for 


Q[p\ix,y) = ksTX^ ■ Xy{j}{x)5{x - y)^. 


(4.22) 


This form will be useful in the following to formally extend the large deviations results that we described in 
Sec.to the present case. 

In the Dean equation, the noise term becomes small when N becomes large. In particular, in the N ^ oo 
limit, one obtains the deterministic non-linear Fokker-Planck equation 


dtp{t, x) + X ■ jp{t, x)=0 (4.23) 

known in the mathematical literature as McKean-Vlasov equation (55]. Some authors refer to the above 
equation as the Vlasov-Fokker-Planck one. In the following, we call the evolution described by this equation 
the mean-field dynamics. 

The above reasoning showed that the McKean-Vlasov equation should describe the evolution of the 
empirical density in the V —> oo limit. This is actually corroborated by a number of rigorous results, since 
the original paper of McKean |66j . More precisely, the property of propagation of chaos was proved under 
mild hypothesis on the smoothness of the b and V, see m and references therein. Moreover, a bound on 
a proper distance between the solution to the McKean-Vlasov equation and the empirical measure at time 
t, depending on the distance at t = 0 is known. Those results are the analogues in the present context of 
the more famous ones due to Braun-Hepp m and Dobrushin |25j for deterministic particles with mean-field 
interactions, see also m- 

The Dean equation suggests that the evolution of the empirical measure for finite but large N is described 
by a weak random perturbation (of order 1/'/N) of the McKean-Vlasov equation. We are then in a similar 
context to that of the Freidlin-Wentzell theory discussed in Sec.[^ except for the fact that the stochastic 
dynamical system is now an infinite dimensional one. We may, nevertheless, hope to obtain large deviation 
estimates in a similar manner by working at the formal level. This will be done in Sec. |4.4| Before, however, 
let us further investigate the McKean-Vlasov equation by considering a simple model system. 
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4.3 Shinomoto-Kuramoto model: the mean-field behavior 


Models in the class of Eq. (4.1) can display a very rich mean-field dynamics. As a simple example, we 


consider in this section the Shinomoto-Kuramoto model introduced in m, describing its mean-field behavior 
by focusing on results that can be obtained analytically or semi-analytically. 

The Shinomoto-Kuramoto model is a one-dimensional model where N particles move on a circle of unit 
radius and are thus identified by their angular coordinate Xn = On defined modulo 27r, for n = 1,..., iV. The 


model is obtained from Eq. (4.1) by setting b{9) = F — h sin(0) and V(9) = (1 — cos 0), where E, h, J are real 


constants that we shall take non-negative. The equations of motion are 


J 


N 


= {F - hs\n9n) - ^ X! Vn , 

m—1 


(4.24) 


where T]n{t) are independent standard scalar white noises. Observe that negative F is related to positive F 
and negative h to positive h by the changes of variables >->■ —and -I- tt, respectively. One of the 

parameters among F, h, J, T is redundant, as it can be fixed by rescaling the other three and time. We could, 
for example, fix the coupling strength J. We prefer, however, to leave all the parameters because we shall be 
interested in limiting cases where one of them vanishes. The Shinomoto-Kuramoto system is related to the 
more famous Kuramoto model |58j of frequency synchronization phenomena in coupled rotators, from which 
it is obtained by setting all the natural frequencies to the same value F and adding white noises acting on 
each rotator, see [1] for a review on Kuramoto and related models. 

For general value of the parameters, the Shinomoto-Kuramoto stochastic dynamics (4.24) breaks the 
detailed balance. This is always true except for E = 0, where the model defines an equilibrium dynamics and 
reduces to a kinetic version of the mean-field ferromagnetic XY model, with describing the angles of planar 
spins and h the external magnetic field h in the X direction. For E > 0, the model may be still interpreted 
within ferromagnetism, except that the planar magnetic field should be taken rotating with angular velocity 
E and the spin angles described in the frame rotating with it. The case /i = 0 is also special since its dynamics 
only trivially breaks the detailed balance that may be restored by returning to the original frame. Thus, the 
Shinomoto-Kuramoto model can be seen as a non-equilibrium version of mean-field ferromagnets. 

for N 


As discussed in Sec. 4.2 


oo, the evolution of the empirical density is described by the McKean- 


Vlasov equation (4.23) which, in the present case, reads: 


dtp + dejp = 0, (4.25) 

jp{t,9) = p{t,9) ^F — hsinO — J J sin{9 — 9) pi^t,!}) — ksT dgp{t,9). (4.26) 

The mean-field behavior of the Shinomoto-Kuramoto model was first studied in m, focusing on long-time 

behavior. In this work, the authors expanded in Fourier modes the stationary McKean-Vlasov equation (4.25) 
and numerically integrated the resulting coupled ordinary differential equations. 

For E = 0, the system relaxes to a stationary solution that is unique for h > 0 and for ft, = 0 and 

ksT > J/2, with p flat in 9 (unmagnetized state) in the latter case. For ft = 0 and ksT < J12 there is a 

one-parameter family of stationary solutions differing by rotation, with p bumped around some value of 9, 
a well known picture for equilibrium ferromagnets of a magnetized state spontaneously breaking the planar 
rotation symmetry. This equilibrium case, also known under the name of Brownian mean field model, is 
studied in detail in miiiiTi]. 

For E > 0, a more complicated phase diagram for the McKean-Vlasov equation emerges, however, see the 
left part of Figure]^ For sufficiently high ft or T, the system relaxes to a stationary solution of Eq. (4.25), 
which for ft = 0 and ksT > J/2 has p flat in 9, as before. For sufficiently low values of ft and T, the long-time 
behavior is, instead, periodic. For ft = 0, the periodic phase sets in for ksT < J/2, as follows from the 
relation to the E = 0, ft = 0 system mentioned before. The periodic and stationary regions are separated 
by bifurcations. The rightmost line (blue dots) is a Hopf bifurcation, while the upper line, a saddle-node 
one. The two bifurcations meet forming a Takens-Bogdanov bifurcation. The careful analysis of bifurcations 
occurring in this model has been performed in m and is confirmed here. In particular, a tiny region where 
the system is bistable is found around the Takens-Bogdanov bifurcation, see the right part of Figure Here 
the McKean-Vlasov equation admits two stable stationary solutions. 

A number of analytical results can be obtained. In particular, we show below that all the stationary 
solutions (stable and unstable) can be studied analytically. Moreover, analyzing their stability we can trace 
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Figure 3: Phase diagram obtained with our semi-analytical results describing the long-time behavior of the Shinomoto- 
Kuramoto model for J = 1 and F = 0.2. We have checked that the phase diagram for other values of the parameters is similar. 
Our results are in agreement with those first obtained numerically in d I73| . Qualitatively, the long-time behavior of the 
Shinomoto-Kuramoto model is very different for high T and/or h and for low T and h. In the first case, the system is stationary 
at long times and the empirical density converges to the unique stationary stable solution of the McKean-Vlasov equation. 
This is true except for a very small region where two stationary stable solutions of the McKean-Vlasov equation are present, 
represented by the shaded region in the inset. For low T and h, no stationary stable solutions of the McKean-Vlasov equation 
exists and the empirical density converges to a periodic solution. The two regions are separated by two lines corresponding 
to the Hopf (dots) and saddle-node (blue line) bifurcations. These two lines merge in a Taken-Bodganov bifurcation. Other 
stationary but unstable states of the McKean-Vlasov equation exist, and will be fully analyzed in the following. 


the bifurcation curves. On the other hand, we were not able to find a closed form for the periodic solutions 
except for the trivial cases where h = 0 or T = 0. 

In Sec. 4.3.1 below, we consider the case where the dynamics respects the detailed balance {F = 0) or 
trivially breaks it {h = 0). Then, in 4.3.2[ we list all the stationary solutions for generic values of the 
parameters and describe how their stability may be analyzed. Finally, we discuss the special case where the 
particles do not interact (J = 0), and the singular situation of zero temperature (T = 0) in Secs. 4.3.3 and 
4.3.4| respectively. 


4.3.1 Equilibrium dynamics (F = 0 or /i = 0) 


We analyze in this paragraph the long-time behavior of the Shinomoto-Kuramoto model for F = 0 or /i = 0. 
We refer the reader to [14] for a more detailed discussion on the case F = 0. 


Let us start with the case F = 0, J > 0 and F > 0. Here, the V-body system (4.24) with finite N respects 


the detailed balance with respect to the invariant measure given by the Boltzmann-Gibbs distribution 


/iv(6»i,...,6»Ar) = ^exp 


1 

kBT 


N 


N 


cos 6»„ - ^ ^ (1 - cos(6'„ - 


(4.27) 


where Z is the canonical partition function. The system is ergodic and the mean of the empirical density 
converges to the one in the Gibbs measure that may be easily calculated. Indeed, applying the the Hubbard- 


Stratonovich transformation to Eq. (4.27), we get the identity 


/7v(^i) • • ■) (^ n ) — 


_ w J 

Ne 

2'kJ ksT Z , 


(m^+h) cose„+m„ — -ij 


dtrij; drUy . 


(4.28) 


The expectation of the empirical density pn{6) in the Gibbs stationary state is equal to the integral of 
/Ar(0, ^ 2 , • • ■, On) over 02 , • • •, that gives 


E 


(2tt)^ ^iVe f w^({m^+h)cose+mySin9-^{Tnl+ml)]i' .yn(h)^\^-^, , //lon'l 

Pn{x) = -- /® ^ ^ duly, duly, (4.29) 
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where m(h) = \/{mx + h)'^ + rriy and 


Z7T 


zi COS 9-\-Z2 sin 9 


dO 


(4.30) 


is the Bessel functions of the first kind. From Eq. (4.29), we can obtain the equilibrium stationary density in 
the N —>■ oo limit by the saddle point calculation: 


PinvW = 1™ E pNiO) = 

N—^OO 


'JT^T d + TUy sin 9 '^ 


e “B 


2ttL 


( m{h) \ 

'{ksT J 


(4.31) 


where {mxjrriy) minimizes 


fir 


,) = 


ml+ml / m(/t) \ 


2JkBT 


‘KkeTJ- 


(4.32) 


The corresponding stationarity equations are 


m{h) 

rrix mx + hriy~kdT 
J mih) 


m{h) 


rr^ 

J 


^ 0 l fc„T I / Iq 


m(h) 


( m.{h) \ 

{ksT J 

( ni{h) \ 
[k^T) 


(4.33) 


with Ii = /q. They may be rewritten as the self-consistency equations 


= J cos 0p2x{d)d9, 


= J / sm9p2x{0)de 


(4.34) 


so that {mx,my) has the interpretation of the magnetization vector (in units of J). 


For h > 0, the stationarity equations (4.33) imply that niy = 0 and nix > 0 solves the equation 


r -\-k, 


kgT ) 
T ! m^+h\ 


(4.35) 


which has a unique solution. In the limit h ^ 0, one recovers the 2"'^-order phase transition (a pitchfork 
bifurcation) located at = J/2. For higher values of T, one has ^ and the magnetization vanishes 

(i.e. nix = O^niy = 0) while for T < T^. the stationary state is spontaneously magnetized {nix > 0,niy = 0). 
Taking h = 0 directly, any rotation of the low temperature solution in 9, with the corresponding rotation 
of the magnetization vector {nix,niy), provides a solution of the saddle point equations which minimizes 


(4.32), whereas nix = 0,niy = 0 corresponding to the flat density gives the maximum of (4.32) instead of the 


minimum. The same limiting densities will be obtained as the stationary solutions of the McKean-Vlasov 
equation for F = 0 , see below. 

Let us now consider the second case, where F,J,T>0 and h = 0. Strictly speaking, the dynamics does 
not respect here the detailed balance. However, it can be recast as an equilibrium dynamics performing the 
change of variables 9{^ = 9n + Ft, which corresponds to sitting on the comoving frame with angular velocity 
—F. The empirical density in this case can thus be obtained from the equilibrium results. Indicating by 
the empirical density for F = 0 and given values of T, J and by p^ the empirical density for the same values 
of T and J but with F 7 ^ 0, we have p^{t, 9) = p‘{^{t, 9 — Ft). We conclude that for T > T^, in the limit of 
long times and large N, the expectation of p^(t, 9) becomes stationary and flat, while for T < Fc it becomes 
periodic. The pitchfork bifurcation at equilibrium is thus modified to a Hopf bifurcation. 


4.3.2 Stationary states for generic parameters 

We now consider generic values of the parameter h,T, J,F with T > 0. To find stationary densities pinv at 


N = 00 , we consider the McKean-Vlasov equation (4.25) in the stationary form, imposing dtPinv = 0. As the 
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Shinomoto-Kuramoto model is one-dimensional, this is equivalent to ask that the mean-field current jp given 
by Eq. (4.26) be a constant, that we shall denote c, for p = pinv We infer that for the stationary densities. 


jpir^^S) = {P - {.rnj; + h) sind -|- ruy cos9]pinv{0) - ksTdepinvid) = c, (4.36) 

where rrix and rUy are given by the equations 

^ 27 r i‘2t: 

rrix = J / cos 9 Pinv (9) d9, rUy = J / sin9 pinv{9) d9, (4.37) 

Jo Jo 

so that {mx,rny) is again the magnetization vector. As the current is constant, jp^^vW — Jpiwi^ + ^o) for 
any 9q implying that 


Jpivv{(^) = ~ [{iT^x + h) cos 00 + Wy sin0o] sin0 -|- [rUy cos9o — [rux + h) sin0o] cos0} pinv{9 + 0o) 

-kBTdgp,nv{9 + 9o) = c. (4.38) 

We can now choose 9q in such a way that 

rUy cos 00 ~ (j^x + h) sin 0o = 0 , {rUx + h) cos 0o -I- rUy sin 0o = y > 0 , (4.39) 


so that Eq. (4.38) becomes 

jpivAd) = {F - y sin0)p„„(0-|-0o) - ksTdgpinviO + 9o) = c. 


Eor any y, the solution to this equation can be written as 


e+2-x 


p ^^^(0 + 0 o ) = j 


dd 


(4.40) 


(4.41) 


with Z the normalization factor. Clearly, for Eq. ( 4.41[) to be a stationary solution, the self-consistency 
conditions (4.37) have to be satisfied together with Eqs. (4.39). Such self-consistency condition can be made 


more explicit by introducing 


Z7T Z7T 

/ h f h 

cos 6 p^nv{0 -\-9o)d9 = J-J cos6»o, fy{y) = / sin9 pinv{99o) d9 = j sin6»o. (4.42) 


0 0 
from which we can eliminate 9q to obtain an equation for y 


[fxiy)- j) +/y(y)^ = ■ 


(4.43) 


We have thus showed that, for any T > 0, stationary solutions of the McKean-Vlasov equation associated 


to the Shinomoto-Kuramoto model (4.25) are given, up to rotation by angle 0o, by Eq. (4.41), where y is 


chosen so that the self-consistency condition (4.43) is satisfied. At the end, angle 0o may be found from 


Eqs. (4.42). The above self-consistency problem can be easily solved numerically. With such procedure, we 
obtained the curves reported in Fig.for the left hand side of Eq. (4.43) as a function of y. Note that such 
a function does not depend on h. Two cases are observed. If T > Ttr, for any value of h only one value of 


y satishes the self-consistency condition (4.43). In this case, whatever the value of h, the McKean-Vlasov 


equation admits only one stationary solution. For T < Ttr^ depending on the value of h, one obtains one or 
three values of y for which the self-consistency is satisfied. Thus, the McKean-Vlasov equation admits one 
or three stationary solutions. Observe moreover that the maximum and the minimum of the curve defines 
two values of y where a pair of stationary solutions is created or destroyed. They thus indicate the values 
of h where saddle-node bifurcations occur. We leave to the reader checking how the self-consistent solutions 
reduce to the ones obtained in the previous section for the equilibrium case F = 0. Let us just note that in 
that case, /^(y) = ^i(^)/-fo(^) and fy{y) = 0. 

To finally recover the phase diagram, we must also analyze the stability of the stationary solutions just 
described. For this, it is enough to look at the linearization of the McKean-Vlasov dynamics around a given 
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Figure 4: We report here the typical curves that are obtained for the left hand side of the self-consistency condition ( |4.43[ l 
that does not depend on h. Two cases are observed. If T > Ttr (on the right), whatever the value of h, only one value of y 
satisfies the self-consistency condition and the McKean-Vlasov equation admits only one stationary solution. For T < Ttr (on 
the left), depending on the value of h, one obtains one or three values of y for which the self-consistency is satisfied. Thus, 
the McKean-Vlasov equation admits one or three stationary solutions. Moreover, the maximum and the minimum of the curve 
defines two values of y for which a pair of stationary solutions is created or destroyed corresponding to the values of h where 
saddle-node bifurcations occur. 


Pinv 5 


dtSp = Rp.^^Sp{t,-), 


Rpi^^Sp = —( 


5p{9) — hsmd — J J sin{9 — d) pinvi'd) dd 

- Pinv{9) J J sin(6* - 9) Sp{9) dd - ksT deSp{9) 


(4.44) 

(4.45) 


One should then study the spectra of the linearized Fokker-Planck operator Rp -^^. This is possible analytically 
only in very special case when pinv does not depend on 9, which holds for h = 0. 

We checked the stability of the stationary solutions by passing to the Fourier transformed picture and 
truncating the resulting system to modes k < K, in which way, Rp^^^ is reduced to a K x K matrix. It is 
then a simple numerical task to find the eigenvalues of such a matrix. Only few modes {K ~ 7) are enough 
to already obtain very accurate results. The stationary states and their stability are summarized in Fig.[^ 
from which we reconstruct the phase diagram anticipated in Fig.|^ 

We conclude by observing that it is possible to calculate analytically where the shown bifurcation lines 
end. Indeed, from the equilibrium solution, we know that the Hopf bifurcation crosses the h = 0 axis in 
T = Jj2. Moreover, it is simple to show that for T = 0 the bifurcation occurs at F — h. 


h 

A: stable 



ao ai Q? <X3 a4 cs 



Figure 5: Stationary states in the mean-field limit of the Shinomoto-Kuramoto model along with their stability. In each region, 
we report which stationary states are present and whether they are stable or unstable, where the letters A, B, C, D indicate from 
which self-consistent solution the stationary state is obtained, see Fig.[^ 
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4.3.3 Free particles (J = 0) 

The stationary solutions for the free particle case ( J = 0) can be obtained very easily directly, 


e+27r 


p_(0) = J 




(4.46) 


or from Eq. (4.411 since in this case mx = niy = 0 = 0 q and thus y = h. 


4.3.4 Zero temperature 

There are two kinds of stationary solutions for T = 0: delta-function-like and smooth solutions. Smooth 
solutions are physically less relevant because they are linearly unstable in the ferromagnetic model (J > 0) 
that we consider here. We discuss smooth solutions in Appendix|^and consider here only the delta-function¬ 
like ones. 

The stationary McKean-Vlasov equation in the form ( |4.40[ ) for T = 0 is solved by 

P^nvi0 + 0o) = S{e - 0i) (4.47) 

if 01 satisfies F — ysm0i = 0. Clearly, such solutions exist only for y > F. The self-consistency condition 


(4.42) reduces in this case to the relation 



(4.48) 


that has no solutions for y if h < F, one solution if h = F and two solutions for h > F. They are 
straightforward to understand. Indeed, one may show that at such solutions, F — /isin(0o -I- 0i) = 0, so 
that the delta-function stationary states describe a situation where all particles sit in the zero of their drift 


F — hsin0, which clearly provides a solution of the equations of motion (4.241 when T = 0. Such solutions 
exist for all N if and only if h > F. For large N, their stability depends on the sign of hcos{0o + 0i), with 
the positive one corresponding to stable solutions. 


4.4 Large deviations for large but finite N 

For large but finite N, the dynamics of the empirical density pN(t,x) deviates from the mean-field evolution 
given by the McKean-Vlasov equation (4.23). Formally, pN solves the Dean equation (4.12) which is an 


infinite-dimensional dynamical system perturbed by a weak noise. Below, starting from this equation and 
employing a functional integral argument, we derive the results obtained rigorously by Dawson and Gartner 
nmn] that describe the dynamical large deviations of the empirical density for mean-field diffusions. We 
also formulate the functional Hamilton-Jacobi equation for the quasi-potential describing the asymptotic form 
of stationary distribution for the empirical densities. Another functional integral argument gives the large 
deviations asymptotics for the distribution of dynamical fluctuations of the empirical current. A generalization 
of the Freidlin-Wentzell theory to infinite dimensional dynamical systems whose deterministic part is given by 
the Fokker-Planck operator with, possibly, non-linear diffusion and drift coefficients is known in the physics 
literature as the Macroscopic Fluctuation Theory. It has attracted much attention in last years in non¬ 
equilibrium statistical mechanics mainly because of its applications to stochastic lattice gases, see [3] for a 
review. Our results may be viewed as an application of that theory to diffusions with mean-field interactions 


of the general form (4.1). It will be also straightforward to specify the results that follow to the case of the 


Shinomoto-Kuramoto model. 


4.4.1 Dynamical large deviations for the empirical density 

Let us fix two times U < tf and consider the probability that the empirical density pNit, x) is arbitrarily close 
to a given trajectory p(t,x) for U < t < tf. In fact we shall be interested only in the large deviations where 
such probability is described by the rate function A[p{-,-)\ as in Eq. (2.7) where trajectories x{t) become 
those in the infinite-dimensional space of densities equipped with an appropriate norm and the limit e —>■ 0 
is replaced by V —> oo. Informally, 


= - lim — lnP[p(-,-)], 

Af-foo JV 


(4.49) 
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where P[p{-, •)] is the probability distribution function(al) describing the distribution of empirical densities 
in the time interval This rate function is a functional of the trajectory p(-,-) and will provide a 

generalization of the Freidlin-Wentzell action of Eq. (2.8) to the present infinite-dimensional setup. 


We may find the form of A[p\ by the Martin-Siggia-Rose functional integral argument [SSj. A similar 
derivation to the one presented below may be found in [53j . Assuming the Dean equation, the distribution 
function P{p\ is that of its random solutions and may be written in the form 


P[p] = E 


5 ( dtp + V • jp -l- 


2kBT 

N 


V • (Tp 0 det 


s(dtp + X^ jp + 4) j 


Sp 


(4.50) 


if we fix the initial value p(ti, x) = Pi{x). The Jacobian factor can be dropped, as it will not contribute to the 
large deviations. Expressing now the delta functional as an oscillatory functional integral in the imaginary 
directions over fields u(t,x) such that / u{t,x)dx = 0, we obtain 


P[P] 


E 


exp 


exp 


N 

2kBT 


N 

2kBT 
‘t 


^ dt 


dtp + V-jp + ,/^V 


N 





■ 

udx 

Vu 


■ jp) u + \p (Vu)^^ dx 


Vu , 


(4.51) 

(4.52) 


where the last equality was obtained by calculating the Gaussian expectation over the white noise ^{t,x). 
The rate function A[p\ can now be extracted by the saddle point argument: 


A[p\ = — mir 


^keT 


{ i^tP + V • jp ) M -f ip (Vm)^) dx , 


(4.53) 


where the minimum is taken over all functions u(t, x) with U < t < tj and spatial mean zero. The minimum 
is easy to calculate resulting in the formula 


ft 


where (—V • p{t, •)V) ^{x, y) is the kernel of the inverse of the operator —V-p{t, •)V in the action on functions 
of zero mean. Eq. (4.53[) may be also rewritten in the form 


tf 


A[p\ = 


(^fdtf {dtp + ^ ■ jp) udx'^ 


AksT Jdt f p{Vu)^dx 

ti 


(4.55) 


which was rigorously derived for the large deviations rate function A[p] by Dawson and Gartner in |181119] . 
Indeed, Eq. (4.55) results from (4.53) upon substituting u ^ Xu and minimizing over A S M. On the other 


hand, expression (4.54) is completely analogous to the Freidlin-Wentzell action (2.8) as it may be rewritten 
as 

Ap] = \ ^^tP - IC[p\){t,x)Q~^[p{t,-)]{x,y) [dtp- IC[p]){t,y)dxdy, (4.56) 


where /C[p] is the drift and Q[p] is the noise covariance of the Dean equation that are given by Eqs. (4.19) 
and (4.22), respectively. 


In our arguments, we ignored the fact that the empirical density pN{t,x) defines a genuine stochastic 
process only after integrating it against a spatial test function. The large deviations for such stochastic 
processes will, however, be governed by rate functions that are contractions of A[p] obtained by minimizing 
it with imposed integrals against the test functions. 


4.4.2 Quasi-potential and Hamilton-Jacobi equation 


The functional Fokker-Planck equation associated to the Dean equation ( |4.12 1 describes the evolution of the 
probability distribution function Pt [p] to observe a density profile p{x) at time t and takes the form 


dtP = Ap, 


(4.57) 
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where £+ is the adjoint of the generator (4.16) calculated with the rule (;^)^ = Analogously to the 

case of finite-dimensional systems, we define the quasi-potential as the functional 


■A[p]=- lim — lnPoo[p] 

A/—>-oo I\l 


(4.58) 


where Poo denotes the stationary solution of the above Fokker-Planck equation. Proceeding formally, all the 
structure of Freidlin-Wentzell theory generalizes to the present case. We are thus very brief in what follows, 
as the derivation of the statements below is very similar to the one for finite dimensional systems, albeit only 
formal here. 

As in the case of finite dimensional systems, we only consider cases where the attractor of the deterministic 
evolution dtp = lC[p\ is a fixed point Pmv The quasi-potential with respect to can be obtained from the 
minimization problem 


T\p] = min A\p ]. 

{p{t,x) I p(-oo,x)=pi„.u{x) , p(0,x)=p(x)} 


(4.59) 


Alternatively, one can try to solve the Hamilton-Jacobi equation 


6A 

Sp{x) 


IC[p]{x) dx 


Q[p\{x,y) ^^dxdy = 0 


5p{x) 


5p{y) 


(4.60) 


that can be obtained by inserting Ansatz Pao[p\ ~ e into the functional Fokker-Planck equation (4.57). 

Equivalently, we can obtain it from the Hamilton-Jacobi equation for finite-dimensional systems (2.24). 


Finally, with the same argument as in Sec. |2.3[ it is simple to show that the minimization in Eq. (4.59) is 
achieved over the solution of the fluctuation or instanton dynamics 


dtp = lCr[p\, K:r[p\{x) = K,[p\{x) + 2 J Q[p]{x,y)dy 


= IC[p]{x) + 2kBT{-Vpix)V) 


SP 

Sp{y) 
5F 


Sp{a 


(4.61) 


4.4.3 Dynamical large deviations for the empirical current 


Empirical current (4.3) is a more general quantity than empirical density (4.2). Indeed, pN and jv satisfy the 


continuity equation (4.4). Given the empirical current Jn in a time interval s-nd the initial condition 

for piy(ti,x) = Pi{x), we can reconstruct the empirical density in the same time interval as pj^(x,t) = 
Pi{x) — /(. V • j]S!{x,s)ds. The converse is however not true, as two empirical currents that differ by a 
divergent-free quantity give rise to the same empirical density. 

Let us consider the probability distribution function P[j(-,-)] of empirical current jN(t,x) in the time 
interval U < t < tf for fixed initial condition for empirical density. The corresponding large deviations rate 
function will be denoted by A’^: 

= - lim ^lnP[j(-,-)] ■ (4-62) 

Af-foo JV 

It may be found using again the Martin-Siggia-Rose functional integral formalism: 


P[j] X E 
= E 


S[j - jp - 


exp 


N 

2 kBT 


dtl (J-Jp- 


J 2^^ J^dtJ (-(j - jp) -A+^pA^^dx 


VA 


VA. 


(4.63) 

(4.64) 

(4.65) 


where A is vector valued, and evaluating the last functional integral with the saddle point method. This 
gives: 

nj-jp? 


1 




dx , 


(4.66) 


where p{t, x) = Pi{x) — /q V • j{x, s) ds. 
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Eq. (4.66) is known in the literature as the fundamental formula of Macroscopic Fluctuation Theory. It 
is often viewed as giving the rate function for the joint probability of observing a trajectory of the empirical 
density and of the empirical current arbitrarily close to p{x,t) and to j{x,t), provided that the latter satisfy 
the continuity equation, and equal to infinity otherwise. This just expresses the fact that density trajectories 
are completely determined once we chose the current trajectory and the initial condition for the density. 

Again, we ignored the fact that these are the integrals of the empirical currents against test functions, 
now both over time and space, that make sense as random variables. The large deviations for such random 
variables will, however, be governed by the rate functions obtained by contractions of that minimize it 
with constraints imposed for the space-time integrals of j. 

Let us conclude by observing that, to the best of our knowledge, large deviations for currents in diffusions 
with mean-field interactions have not been discussed in the mathematical literature. We will devote a future 
publication to the investigation of current fluctuations in the Shinomoto-Kuramoto model. 


4.5 Perturbative calculation of the quasi-potential 


In this final section, we discuss the perturbative calculation of the quasi-potential for the diffusions with 
mean-field interaction. In Sec.[^ we developed a perturbative scheme to calculate the quasi-potential for 
finite dimensional systems. The first objective is then to translate that perturbative scheme to the infinite¬ 


dimensional setting described by the Dean equation (4.18), which we do in Sec. 4.5.1 


It is possible to obtain explicit results in two cases: close to the free-particle dynamics J = 0, see 
Sec. |4.5.2] and for the Taylor expansion of the quasi-potential around stationary solutions of the McKean- 
Vlasov equation, see Sec. 4.5.5] We shall see that in these cases the perturbative expansion of the quasi¬ 
potential reduces to solving partial differential equations instead of functional differential equations. 

For the Shinomoto-Kuramoto model, in the case of perturbations close to the free-particle dynamics 
J = 0, we present explicit results at the I®* order in J, see Sec. 4.5.3 Interestingly, the quasi-potential that is 
a local functional of p for the unperturbed J = 0 case, is shown to become non-local already at the I®*^ order 
in J. The resulting analytical expression is evaluated with a simple numerical scheme and we also compute 
the rate function for the fluctuations of single-particle observables such as the magnetization. Moreover, we 
discuss how the numerical scheme can be generalized to higher orders stressing that this does not increase 
the computational complexity. 

We also discuss the Taylor expansion of the quasi-potential around the fixed points of the McKean-Vlasov 
dynamics. It would be possible to obtain explicit results also in this case. We detail how an algorithm can 
be designed for this purpose in Appendix]^ We did not, however, perform the corresponding calculations, 
leaving them to future investigations. 

Although the perturbative techniques developed in Sec. 3 in the context of finite-dimensional systems are 
easily transposable to the infinite-dimensional setup discussed here, we are unable to say anything about the 
nature of the resulting perturbative expansions (are they convergent? asymptotic?) as, to the best of our 
knowledges, the techniques used in finite dimensions to settle such questions do not extend in a Banach space 
setup. 


4.5.1 Perturbative schemes in infinite dimensional setting 

We generalize here the perturbative schemes for calculating the quasi-potential, developed in Sec.[^ to the 
case of infinite dimensional systems at hand. This can be done very easily either by starting again from the 


Hamilton-Jacobi equation, now in the infinite dimensional version (4.60), or by formally taking the limit of 


infinite dimensions in the expressions of Sec.[^ Such a formal limit transforms derivatives into functional 
derivatives and finite-dimensional scalar products into ones. 

Both the expansion centered on the attractor p^^^ of the perturbed system (a stable fixed point of the 
McKean-Vlasov dynamics) and the direct one can be easily obtained. We only consider here the second one, 
as it will be used in the following. We denote the direct expansion of the quasi-potential as 




(4.67) 


dropping the hats of Sec. 3.2 


and analogously for JC^[p] and Q^[p]. 
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With the procedures outlined above, the hierarchy (3.21) and (3.22) is now replaced in the infinite dimen¬ 
sional setting by 


5p{x) 


r 


dx = 0 


where ici^^ is the drift of the fluctuation dynamics for the unperturbed problem: 


0 

for n = 0 , 

(4.68) 


for n ^ 0, 

(4.69) 


4°)[p]=/C(°)[p]+2Q(°)[p]- 


Sp 


(4.70) 


n—1 


and is given by 

k^l 

n—k 

+ E 

1=0 


SJ7(n-k) 

5p{x) 

§jr{n-k-l) 

5p{x) 




Q^'"Hp]{x,y) dy 


dx 


dpiy) 


6^ 

Sp{x) 


r 

^(-)[p](x)+ J QM[p]{x,y)j^^dy 


dx. 


(4.71) 


As it will be clear, there is a full analogy with finite-dimensional systems and we refer the reader to Sec.[^ 
for the discussion on the perturbative schemes. In particular, we recall that the solution of these equations 
obtained with the methods of characteristics is unique and an explicit expression may be written in terms of 
the solution for the 0*^*'-order fluctuation dynamics. In the actual infinite dimensional setting, this solution 
is given by 


7° 

J —oo 


where p(t, x) is the solution of the fluctuation dynamics 

dtp = /C(°) [p] 


(4.72) 


(4.73) 


starting at t = —oo at the stable fixed point p^^^ = pmv of the unperturbed McKean-Vlasov dynamics and 


satisfying the final condition p(0,a;) = p[x), compare to Eq. (3.27). It should be clear from the discussion 


of Sec.[^ that the convergence of the above integral in this infinite-dimensional setting is assured once the 
unperturbed fluctuation dynamics (4.73) linearized around pmv escapes fast enough from pmv This happens, 


for example, in the typical case in which its generator has a spectral gap. We finally recall from Sec. |3.2| that 
the constants may be iteratively fixed but do not enter the expression for 5*^"^ ... ^jr(n-i)] gQ 

largely irrelevant and may be adjusted conveniently. 

4.5.2 Expansion around the free particles dynamics (J = 0) 

We now consider the perturbative expansion of the quasi-potential T around the independent particle dy¬ 
namics, choosing the perturbative parameter A = J. In this case, 


^■'^(a:) = -V 


p{x){h{x)-J J{W){x - y) p{y) dy'^ - ksTVpia 


(4.74) 


and 


Q'^[p](.x,y) = ■ Vy{p{x)S{x - y)) . 


(4.75) 


For J = 0, particles do not interact. Even if for J = 0 the system breaks the detailed balance, the quasi¬ 
potential is known here as its explicit form follows from Sanov’s theorem: 


= I 


/ \ 1 / \ / \ 

p{x) In-— - p{x) -f Pinv[X) 

PinvKp^) 


dx , 


(4.76) 
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where pmv is the stable stationary solution of McKean-Vlasov’s equation which reduces here to the linear 
Fokker-Planck equation for a single-particle diffusion. Hence pinv exists and is unique under simple assump¬ 
tions about the single-particle drift b(x). In the case of the Shinomoto-Kuramoto model, we have explicitly 
calculated pinv for J = 0, see Eq. (4.46). We observe that, despite the fact that the dynamics breaks in 
general the detailed balance, the quasi-potential for J = 0 is a local functional of p. 

As anticipated, it is convenient to perform the direct perturbative expansion, for which we have 


and 


/C(")[p](cc) = 

Q^°'>[p\{x,y) 

Q^'^\p]{x,y) 


V • [-p{x) b{x) + ksTVy^x )], 


V 

0 


p{x) J {VV){x-y)p{y)dy 
for n >2 

ksTVx ■ Vy{p{x) S(x - y)) 
0 for n > 1. 


(4.77) 

(4.78) 

(4.79) 

(4.80) 

(4.81) 


It should be clear now why, from a practical point of view, it is simpler to deal here with the direct expansion 
than with the expansion centered on the attractor of the perturbed dynamics. Indeed, in the former case, 
almost all and 2^"^ are zero. This would not be true in the attractor-centered expansion because the 
attractor of the perturbed dynamics depends on J. 


The quasi-potential at the 0*^ order is given by Eq. (4.76). We thus have 

= In ' ' 

,(a;) 


_=ln^ 

(5p(ir) Pinvi 


The 0*^ order fluctuation dynamics takes the form dtp = /Cr°^[p], where 

[p] = - V • (p 6) + 2kBT V • [p V In p„„] - V^p = K^o^p. 


(4.82) 


(4.83) 


Note that is a linear operator. When detailed balance is respected (i.e. b = —VU for some potential 
U), one can easily check that 


vW[p] = -V • (pVV) - ksTV^p = -/c(°)[p] 


(4.84) 


where we have used that pmv = {1/Z) exp{—U/ksT) in this case. The above relation just reflects the fact 
that the fluctuation dynamics is the time-reversal of the relaxation dynamics for equilibrium problems. 

Let us specify the hierarchy ( |4.69 ) to th e pres ent situation. The left hand side is clearly unchanged, 
except that k!^\p] takes the particular form (4.84), while the right hand side becomes equal to 




(4.85) 


for n = 1, and to 




X-pin-k) AF(fe) 

(4,86) 


for n > 1. A very important point is that, in this case, the solutions of the hierarchy are for n > I, 

homogeneous polynomials of degree (n -I-1) (up to constants), 


J'l"') [p] = -f -— {xo,...,Xn) p{xo) ■ ■ ■ p(x„) dxo ■■■ dxn , 

(»^ + 1)! J 


(4.87) 


where the kernels are symmetric in {n + 1) variables and the choice of constants is essentially 
irrelevant. Using Eq. (4.87), we infer that the left hand side of the hierarchy (4.69) can be written as 


5p{x) 


lC^^\p]{x) dx = ^ J (I>^'^\xo, . . ■,Xn)J p(xo) ■ ■ ■ p(Xn) dXQ ... dXr, 
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(n+ 1) 




In the above expressions, denotes, in a slightly abusive notation, the adjoint operator to K^\ 

= {h- 2kBT (V lnp„,)) • V, - fc^TV^ , 


see Eq. (|4.83|), in the action on a function of variable xq- Similarly, n+i denotes the operator 


(4.89) 


^(0)t _ /r(0)t 

r,n+l / ^ ^^r,Xm 


(4.90) 


m—0 


acting on functions of n-tuples (xq, ..., a;„). Note that the differential operator Kr^'^\x) is the generator 
of a diffusion process with drift —b{x) + 2kBT{W In pmv) which is the time-reversal of the original single 
particle diffusion. This operator has a one-dimensional kernel composed of constants and a one-dimensional 
cokernel composed of functions proportional to Pmv For the same reasons, the operator Fas also a one 

dimensional kernel composed of constants and a one dimensional cokernel composed of functions proportional 
fO Pinv (xo) • • • Pinv(x„). We shall use those properties below. 

In order to compute the right hand side of the hierarchy, we shall treat separately the cases with n = 1 
and with n > I. For n = 1, we have 

5(l)[jF(0)] 

= - J (V^E)(x - y)-I-|(VE)(x - y) • (Vlnp„„(x) - Vlnpi„„(y)) p{x) p{y) dx dy , (4.91) 


where we have made the integrand symmetric under the exchange of x and y. Comparing Eqs. (4.881 and 
(4.91), we obtain a differential equation for 0^)^ 


(^r° 2 V^^^)(a:>y) = -2(V^E)(x - y) - (Vy)(x - y) • (Vlnp™„(x) - Vlnpi„^(y)) . 


(4.92) 


For n > 1, the right hand side of the hierarchy can be written in terms of the kernels (j)^^\xo, ■. ■ ,Xk) as 
^ (^3:o^)(^0 ^n) ’ ^)(^05 ■ ■ ■ 5 ^n—l) piX{)) ‘ ’ ’ piXn) dxp • * * dx^i 

Tl 1 J n 

~ X] k'.{n- k)\ J ■■■,Xk) ■ (Vx^<()^”“''^)(xfc,... ,x„)p(xo) • • • p{xn) dxo ■ ■ ■ dxn ■ (4.93) 

Then, for any n > 1, we obtain the differential equations 


(^ 0 , ..., X„) = (n -k 1)! Sym 


h {^xoV)ixo - x„) • ^)(xo,. . .,X„_i) 


n—1 


- ksT E ^^;^(V,,.^('=))(xo,..., Xfc) • (V,,(^("-'=))(xfc,..., x„) 


k=l 


(4.94) 


wher e Sym [ • ] stands for the symmetrization of the argument with respect to variables xq, .. •, x„. Eqs. (4.92 1 
and (4.94) may be solved iteratively for the kernels ^^"■^(xo,..., x„) provided that their right hand sides 


are orthogonal to the one-dimensional cokernels of operators spanned by the product functions 

Pinvixo) ■ ■ ■ pinv{xn)- Thesc are the solvability conditions S^'^'>[pinv\ = 0 that were discussed in Sec. 


3.2 


in the finite-dimensional context, see Eq. (3.24). Such conditions are easy to check in the lower orders but are 


difficult to prove order by order. They are, ho weve r, implied by the existence of the perturbative expansion 
centered at the attractors, as discussed in Sec. 3.2 Observe finally that are defined by the above equa¬ 


tions only up to constants that are in the kernel of This ambiguity just leads to shifts in constants 


C(0 

appearing in Eq. ( |4.87[ ) which are not relevant (remember that densities p are normalized). 

We have greatly simplified the problem of calculating the quasi-potential perturbatively around the free 


particles dynamics. Instead of solving the functional diff erent ial equations (4.69) we have to sol ve the partial 
differential equations for the kernels defined by Eq. (4.87). These differential equations are (4.92) for the 
fst Qj.(jg]- contribution and (4.94) for the higher orders. Solving them is a rather simple numerical problem. 


In the next section, we discuss the results obtained in the case of the Shinomoto-Kuramoto model. 
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4.5.3 Expansion around J = 0 : results for the Shinomoto-Kuramoto model 

Let us now specify the discussion to the case o f the Shinomoto-Kuramoto model introduced in Sec. |4.3| and 
describe some explicit results. In this case, Eq. (4.92) for the reduces to the relation 


where 


and 


=-2cos{9 - §) - sm{9 - ■&) ■ [dglnpinviO) - d^lnpinvi'^)] , 
X(o)t = [F-hsm9- 2kBT {de In p„„)] dg - keTdl . 


(4.95) 


(4.96) 


(4.97) 


We observe that identity (|4.95 ) is a Lyapunov equation and several techniques can be employed to solve it, 


see for example [56l 0 I69| l The one we have chosen is to perform an expansion on the eigenfunctions of 
operators and 


Let us first set the notations. For k = 1,..., we denote by Uk the eigenfunctions of and by Vk 
those of As is not self-adjoint with respect to the scalar product, these eigenfunctions are not 

connected by complex conjugation but we may assume that the corresponding eigenvalues are: 




'^k — ^k'^k 


and 
do) 


Vk = ak Vk 


(4.98) 


where oik denotes the complex conjugate of ak- Ki- is a Fokker-Planck operator. Zero is its simple eigenvalue 
and pinv is the corresponding eigenfunction. We choose it to be equal to ui- Other eigenvalues of Kr^'^ have 
strictly positive real parts. Similarly, we choose Ui = 1. It corresponds to the unique zero mode of 
We can also assume that the eigenfunctions are mutually orthonormaj^ in the following sense 


Uk{0) vi{9) d9 = 6ki ■ 


(4.99) 


It is then simple to show that 


</.«( 0 ,^?) = 


5. 


( 1 ) 

kl 


V 


Uk{O)ui{'0) , 


(4.100) 


where 




(4.I0I) 


(note that a?; -I- a; ^ 0 for (fc, 1) ^ (I, I)). One has S'!! = 0 and we have chosen to be orth ogona l to the 
kernel of composed of the constant functions. The calculation of constant of Eq. (4.87) will be 


discussed afterwards. 


Eq. (4.100) permits to obtain numerically the quasi-potential at the 1®*^ order in J. One first looks for 


an approximation of the eigenfunctions and eigenvalues of and of by expanding Eqs. (4.98) in the 


Fourier modes and truncating the hierarchy in order to deal with matrices. The calculation of through 
Eq. (4.100) is then easily performed. We checked that very few Fourier modes have to be retained (> 20) to 


obtain an excellent approximation of eigenfunctions and eigenvalues. Even less eigenfunctions (> 5) need to 
be employed in the calculation in order to obtain quite accurate results. 

Although the scheme proposed here is a bit more involved than other methods, such as a direct Fourier 


expansion of Eq. (4.95), it is much more powerful than the latter. Indeed, in the method proposed here, we 


only need to diagonalize the matrices that approximate operators Kr^'^ and that act on the space of 


functions of one periodic variable. Instead, attempting to directly solve Eq. (4.95) by expanding on Fourier 


modes results in diagonalizing the operator acting on the space of functions of two variables. Even 

more important is the fact that the method proposed here is fully generalizable, without increasing the com¬ 
putational complexity to obtain higher perturbative orders (encoded in kernels for the quasi-potential. 


^This is possible whenever the eigenvalues are all distinct. We have checked numerically that this is indeed the case. 
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Figure 6 : On the left: the kernel kBT4>^'^\9,'d) that permits to obtain the first-order correction to the quasi¬ 
potential of the Shinomoto-Kuramoto model in the perturbative expansion around the free particle dynamics 
with J — {). The values of parameters are F = 0.2, T = 0.3, h = 0.5. This figure was obtained by numerically 
solving Eq. (4.95) as detailed in the text. On the right: the same function but in a contour plot. 


In contrast, in the direct Fourier expansion, one needs to diagonalize matrices of dimension {kmax)'^ where 
kmax is the number of retained Fourier modes, making the problem practically intractable already for small 

n. 

We report the results for the 1®* order correction to the quasi-potential in Fig.|^ for some typical 
choice of parameters for which the McKean-Vlasov equation has a single stationary stable solution (namely 
F = 0.2, T = 0.3, h = 0.5). A rather non-trivial result emerges. As a check of the accuracy of the 
perturbative expansion, we also compared the exact stationary state of the system for J > 0 to the one, 
Pinv,i, obtained by imposing that the functional derivative of 


J-(O) + jjr(l) ^ 


p{e) de - 1 


(4.102) 


vanishes. This condition implies, with the use of (4.82) and (4.87), that 


Pinv,l{P^ — ^ 


-J J P^nv,lW 


(4.103) 


d 

inv 


where Z = e'^ fixes the normalization of pinv,i- Recall that the true mean-field stat ionary density p'l. 
minimizes the quasi-potential F"'[p\ with the value at the minimum equal to zero. Eq. (4.103) may be very 
easily solved numerically with an iterative scheme, once is known; we have verified that this procedure 
converges after few iterations. In Fig.[7) we also compare Pinv,i with the exact result obtained from Eq. (4.41), 
which shows that already at the 1®* order the stationary state for the system with J > 0 is m uch b etter 
approximated by pinv,i than by the stationary state at J = 0. As for the constant from Eq. (4.87) with 
n = 1 , it will be convenient to fix it so that the quasi-potential calculated to the 1 ®‘ order -F hag 

the minimal value equal to zero, as does the complete quasi-potential. This is achieved by setting 




(4.104) 


where on the right hand side is taken homogeneous in p (i.e. with set to zero). 

Once the quasi-potential at order J is known, we can also compute the large deviations rate functions 
for the one-particle observables ^ where 5 is a given function. In the following, we consider 

the particularly relevant examples of the magnetization along x and y axis corresponding, respectively, to 
g = cos and g = sin. What follows may be straightforwardly generalized to any choice of g. We denote by 
Xx{(j) (resp. Iy{a)) the rate function for the probability of observing a value of magnetization rux ~ cr (resp. 
my ~ a). Let us consider (the case of Ty may be analyzed similarly). To the 1®* order in J, we have 


IxiA = min 


(4.105) 


where the minimum is taken over all positive functions p with unit integral respecting the constraint 
J cos(d) p{6) dO = fj. The profile realizing the minimum will be called popt in the following. 
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Figure 7: Stationary state of the Shinomoto-Kuramoto model for J = 0.3 obtained from the perturbative 
expansion of the quasi-potential (blue continuous curve) compared to the exact result (red dashed dotted 
curve) given by the Eq. (4.41) and to the stationary state for the system with J = 0 (green dashed curve). 
The values of parameters are F = 0.2, T = 0.3, h = 0.5 on the left and F = 0.2, T = 0.2, h = 0.15 on the 
right. The plots show that the 1®* order correction to the stationary state strongly improves the prediction. 


The minimization problem in Eq. (4.105) can be solved introducing two Lagrange multipliers 71 and 72 


associated, respectively, to the total mass and to the magnetization constraints. A simple calculation shows 
that the minimizing profile satisfies the relation 


Popt(^) — ^ Pi'. 


,{e) exp 


-72 cos 6 »-J j popt{i))d'd 


(4.106) 


where Z normalizes Popt- As for the case of the most probable state, the above self-consistent equation can 
be solved iteratively. We first fix 72 and iteratively solve Eq. (4.106) normalizing the iterated solution to 


unity at each step. Once the optimal profile has been obtained, we then calculate the corresponding value 
of rux- As the initial condition we used the optimal profile for J = 0. We have checked that the iterative 
scheme converges in few steps. The results obtained for and Xy are plotted in Fig. ||for the same choice of 
parameters as those employed in Fig.[^ In the plots, the results obtained at the 1®‘ order in J are compared 
with those obtained for J = 0 in the whole range of possible magnetizations. We observe that the 1®* order 
corrections are significant. Fox example, the probability of fluctuations of rrix smaller than the typical value 
is diminished with respect to the J = 0 case but the probability of very rare fluctuations is instead increased. 


We conclude by noting that, as already discussed, it would be simple to extend the calculations of the 
quasi-potential and of rate functions of observables to higher orders by computing with the method employed 
in this section the higher kernels We do not pursue, however, this direction here, as our scope was 

rather to show that the perturbative approach can lead to explicit results than to analyze in details the 
Shinomoto-Kuramoto model. 


4.5.4 Expansion around the equilibrium dynamics {F = 0 or h = 0) 


When the drift b{x) may be written as b = —VI/ for some U, the quasi-potential can be given explicitly 
because detailed balance holds. We note that the Dean equation (4.18) corresponding to this case has also 


been derived in [HHS], where it is called stochastic Smoluchowski equation. One has in this case 

•^[p] = + ^ / y (x - y) p{y) dy^ dx + J p{x)\n p{x) dx . (4.107) 


We did not develop the perturbative expansion close to the detailed balance because obtaining explicit results 
would be quite complex. Indeed, the equations for the quasi-potential at order n are functional differential 
equations and not, as in the case of a perturbation close to J = 0, ordinary differential equations. Even though 
results may be explicitly obtained by a discretization, the problem is not computationally straightforward. 


4.5.5 Taylor expansion close to a stationary state 

Let us discuss now the Taylor expansion of the quasi-potential F around a stable stationary solution of 
the McKean-Vlasov equation. As we have seen in Sec. 3.3 for finite dimensional systems, this expansion may 
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Figure 8 : Rate functions for the probability of fluctuations of the x (on the left) and y (on the right) 
component of magnetization. The results of our perturbative analysis to the 1®* order in J (blue continuous 
lines) are compared with those of the unperturbed model with J = 0 (green dashed lines). The value of 
parameters used are the same as in Fig.|^in the upper plots while they are set to F = 0.2, T = 0.2, h = 0.15 
in the lower plots. 


be viewed as a particular case of the perturbative expansion of the quasi-potential for a parameter-dependent 
system. In this setting, the role of the unperturbed quasi-potential is played by the quadratic approximation 
to F that, consistently, will be denoted in the following by Then, once j-(o) 

is known, any higher 


order correction may, in principle, be calculated. After extending here the discussion of Sec. 3.3 to a general 
diffusion with mean-held interaction (4.1), we discuss in Appendixhow explicit results could be obtained 
in a similar way as for the perturbative expansion around J = 0. 

We would like to calculate the quasi-potential F[p] for p = {pmv -I- p) as a power series in p. Note p has 
then to have the vanishing integral. In order to use the same notations as in the hnite dimensional setting, 


see Sec. 3.3 we introduce the A-dependent system dehned by 


lC^[p] — —K. [pinv + Ap] 

[p] = Q [Pinv + Ap] 

[p] — \_Pinv Ap] . 

For any n we thus have 

[p] = -V • [p 6 - Jp,nv (VR * p) - jp (VR * P,™)] + keTW^P , 
/CW[p]= JV-[p(VI/*p)], 

A*'"^[p] =0 for n > 2 . 

Clearly, /C^°^[p] = Rpir^^p, where denotes the linearized Fokker-Planck operator. Moreover, 

g( 0 )[p] = -fcBTV-p,„,V, 

Q^^'>[p] = -kBTW ■ pW, 

[p] = 0 for n>2. 


(4.108) 

(4.109) 

(4.110) 


(4.111) 

(4.112) 

(4.113) 

(4.114) 

(4.115) 

(4.116) 


Recall that the computation of the quadratic order of F corresponds, in the above A-dependent system, to 
the calculation of Moreover, with these notations the hierarchy of equations to so lve in orde r to obtain 
the Taylor expansion of the quasi-potential close to pmv is given precisely by Eqs. (4.68) and (4.69) with 
and expressed by the above formulae. 


We shall search for the solution of this hierarchy in the form of homogeneous polynomials 


•^^"^[p] = (^^2)! / ‘P"’i^O,---,Xn+l)p{xo) ■■■ p{Xn+l)dxo 


dx- 


n+1 5 


(4.117) 
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where ip'^{xo,... ,Xn+i) are symmetric kernels on which we shall impose the relations 

J <p"(a;o,..., Xn+i) dxo = 0. 


(4.118) 


As in the other cases, once the quadratic o rder is known, the higher orders can be obtained by using 
the general solution to the hierarchy, see Eq. (4.72). We thus concentrate here only on We introduce 

the operator $ 

•^[pKx) = J (p°{x,y) p{y)dy (4.119) 

that will be considered as acting in the space Hq of functions with vanishing integral and scalar product. 
Note that 4) is a symmetric operator. From Eq. (4.68), we obtain the operator identity 

$i?p_ = 2fcBT$(V • p™„v)$, (4.120) 

where Rinv is also viewed as an operator in Hq. Assuming that 4> is invertible, as it must be from the fact 
that is positive definite, we infer that 




-1 




= 2fcBT(V • 

Pinv V). 


which is the infinite dimensional analogue of Eq. (2.38). Its solution is given by the relation 

poo 

'-1 = -2kBT / e*^" 

^0 




(V • 


dt , 


(4.121) 


(4.122) 


see Eq. (2.37). Similarly to the finite-dimensional case, convergence of the integral is assured if pinv a non¬ 
degenerate stable stationary solution of the McKean-Vlasov equation, i.e. if the spectrum of on the 

space of functions with vanishing integral in contained in the complex half-plane with negative real part. 
Once we know the invertible operator 4?“^, the function (xq,xi) is extracted as the kernel of its inverse. 
With an argument analogous to the one used in Sec.|3.4|for the finite dimensional systems, one can show 


that Eq. (4.122) implies that the covariance of density fluctuations diverges with the mean-field exponent.when 
approaching a codimension-one bifurcations of the McKean-Vlasov dynamics. We note that the latter result 
was obtained rigorously in [5]. 

Let us now consider the 0*^-order fluctuation dynamics. We have 


2 fcBTV • 


(4.123) 


is again an operator in Hq that is invertible and has sp ectrum in the half-pla ne with positive imaginary 
part (it should not be confused with considered in Sec. 4.5.2). With Ansatz (4.117), the left hand side 
of Eqs. (4.69) takes now the form 


5 j:(n) 

5p{x) 


[p] {x) dx 


{n + 2) 


ij pixQ) ■■■ piXn+l)dXQ---dXn+l (4.124) 


in the notation of ( |4.90 ). As for their right hand side, it follows from Eq. (4.71) that 

^7^ / ^"'i^°^---^^''~^^(xQ,...,Xn+l)p{Xo) ■■■ p{Xn+l)dXQ---dXn+l, 


{tl -)- . 


(4.125) 


where the symmetric kernels ^](a:o,.. • ,a;„+i) are expressed in terms of with k < n and 

again satisfy the constraint (|4.118 ). Thus in terms of the kernels, Eqs. (4.69) become the identities 


(■f^r°i +2 <7’") (a^o, • ■ •, a;„+i) = ..., ^](a;o, ..., x „+ i ){ xq , x„+i) 


(4.126) 


which may be solved iteratively since the operators are invertible on (note that the con¬ 

straint (4.118) characterizes the symmetric functions in 77®^"^^^). Further details are left to the reader, 
but in Appendix]^ we briefly discuss how as well as higher order kernels, could be explicitly calculated 
numerically, leaving out of this work the explicit numerical implementation of the proposed algorithm. 
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5 Conclusions 


The main aim of this paper was to discuss a perturbative approach to the calculation of the quasi-potential 
in parameter dependent stochastic dynamical systems. As we have seen, the Taylor expansion of the quasi¬ 
potential around an attractor of the deterministic dynamics was also covered by our theory. 

The approach developed in the paper may be summarized as it was already done in the introduction. The 
perturbative expansion breaks the loop connecting the quasi-potential and the instanton dynamics, with one 
needed to obtain the other. At any perturbative order the quasi-potential can be computed just from the 
knowledge of the instanton dynamics for the unperturbed problem and the lower order results. Indeed, our 
approach gives explicit formulae that permit to iteratively calculate any order of the power series expansion 
of the quasi-potential in the perturbative parameter. Explicit results can be obtained once the instanton 
dynamics for the unperturbed problem may be analytically or numerically solved. 

We have first developed the theory for finite dimensional systems, where the mathematical details can be 
easily handled precisely. Our strategy was to perturbatively solve the Hamilton-Jacobi equation which, as it 
is well known [32) . has the quasi-potential as the unique non-trivial solution under the hypothesis that the 
later is smooth enough. It is also known m that, at least in a neighborhood of the attractor, this is indeed 
the case, while further away singularities may occur, see for example (52 HU EH [53 EH EZ] ■ An investigation 
on how to deal practically with singularities of the quasi-potential is left for the future. We also note that the 
perturbative study of transition rates between basin of attraction requires a specific approach, as was done 
for instance in the reference using Melnikov’s method for perturbations of Hamiltonian systems. 

In the second part of the paper, we moved our attention to many body systems. In particular, we 
considered N particles that undergo an overdamped diffusion, interact with mean field conservative forces, 
and are driven out of equilibrium by an external drift, see Eq. (4.1). In this case it was possible to (formally) 
derive a fluctuating hydrodynamics, with the noise term proportional to I/'/N, that describes in a closed way 
the evolution of the empirical density of particles. This fluctuating hydrodynamics is known in the literature 
under the name of the Dean equation |22j . In the limit when N ^ oo, the empirical density obeys then a 
deterministic PDE, the McKean-Vlasov equation [66] • This provided a formal but quick way to recover the 
results on propagation of chaos that appeared long ago in the mathematical literature. 

In the case of a specific mean-field system, the Shinomoto-Kuramoto model of coupled rotators, we 
showed that such an approach can be employed to obtain explicit results. We first discussed MacKean- 
Vlasov’s dynamics describing the behavior of the system in the iV —>■ oo limit, characterizing analytically or 
semi-analytically all the stationary solutions of the system and their stability. Already this very simple out- 
of-equilibrium mean-field system displays quite a complex phase diagram with bifurcation lines and several 
attractors (stationary states as well as limit cycles). Albeit the Shinomoto-Kuramoto model has been studied 
in literature by numerical techniques, our analytical results characterizing all its stationary states are original. 

After discussing the N = oo limit, we concentrated on large deviations around it. This was done by 
studying the Dean equation for N large but finite. Applying formal techniques from field theory (the Martin- 
Siggia-Rose representation and the WKB approximation), we showed that fluctuations of the empirical density 
are described by a generalization of the Freidlin-Wentzell theory to infinite dimensional systems. This per¬ 
mitted to recover formally rigorous results that already appeared in mathematical literature [mill! and to 
make a connection with the Macroscopic Fluctuation Theory [6]. 

We subsequently generalized our perturbative approach to the calculation of quasi-potential to the case 
of fluctuating mean-field hydrodynamics. For the Shinomoto-Kuramoto model, we explicitly computed the 
quasi-potential close to the uncoupled particle dynamics to the 1®* order in the mean-field coupling strength 
J and found the corresponding rate function for the probability of fluctuations of single-particle observables 
such as the magnetization. Our results seem to provide the first explicit calculation of the quasi-potential for 
diffusions with mean field interaction driven out of equilibrium. Higher orders in J could be also computed 
without much difficulty. As our aim here was mainly to develop the general theory and to illustrate that it 
can be employed to obtain explicit results, we have not pursed such calculations further but discussed instead 
how one could similarly obtain explicit results for the Taylor expansion of the quasi-potential around the 
stationary states. 

The perturbative approach developed in the present paper seems quite general and we expect it to be 
useful in several other problems. For example, it may be employed in the future to study other long-range 
interacting systems which are driven out of equilibrium by different mechanisms, such as stochastic forces 
which are spatially correlated. This setting is natural for systems whose constituents exhibit long-range 
interactions (such as plasmas or gravitational systems) but are exposed to stochastic external fields (electric 
for example) mainly acting on a given spatial scale. The typical behavior of this class of system was analyzed 
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with kinetic theory in |68L I69j and numerical simulations have shown how the phase transitions are modified 
when departing from the detailed balance |t)9j . The perturbative approach exposed here may give an analytical 
insight to such problems. This direction is presently pursued by one of the authors. 

Similarly, applications of such perturbative methods could be found in the study of large deviations for 2d 
turbulence problems, which are relevant for climate modelling liiiniiziiHiiss]- Even though these are quite 
academic questions [S] because realistic models of 2d turbulence and those arising from climate models are 
far from the perturbative regime, this approach can give valuable physical insight as it is difficult to handle 
the rare events in realistic systems. 

Finally, while large deviations for the empirical density of mean-field interacting diffusions appeared in 
literature long ago, mathematical literature does not cover fluctuations of the empirical current. The formal 
approach through the Dean equation presented in this paper clearly permits to go in this direction. Some 
results on the current fluctuations that we have obtained already, as well as further explicit results on the 
calculation of the quasi-potential for the Shinomoto-Kuramoto model and their comparison with results 
obtained from direct numerical simulations are left for a future publication. 
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Smooth stationary solutions of the Shinomoto-Kuramoto model 

at T = 0 


In Sec. 4.3.4[ we have discussed simple delta-like stationary solutions of the McKean-Vlasov dynamics for 
the Shinomoto-Kuramoto model at T = 0. They described the particles accumulated at the fixed point of 
single particle drift. These are however not the only stationary solutions for vanishing temperature and also 
smooth solutions are present, which, however, are linearly unstable in the ferromagnetic model considered in 
this paper (J > 0). We briefly discuss such smooth stationary solutions at T = 0 and their stability in this 
appendix. 

For r = 0, Eq. (4.36) implies that 


PinviP T ^o) — 


27r (l — sin 0) 


provided that 0 < y < F and the self-consistency relation (|4.43 ) holds taking now the form 


& 




(A.l) 


(A.2) 


The left hand side of the last equation grows monotonically from 0 at y = 0 to (F/J)'^ + 1 aX y = F. Hence 
there is a unique solution for 0 < y < F if and only if F'^ + J'^ > h'^. We infer that there is a unique stationary 
smooth solution at T = 0 when the latter constraint holds and none otherwise. 

To examine the stability of such solutions, we examine the eigenfunctions Sp^{0 + 9 q) with vanishing 
integral of the linearized Fokker-Planck operator of (4.44) at T = 0. We may write 5p^{9 + 9 q) = 


deSffj,{9), where Sf^ is periodic on [0,27r]. The eigenequation takes the form 

{Rpi„^deSff,){9) = - ysin6») de5f^{9) - Yi pinvi.6 + 9 q) sinS -f Y 2 Pinv{d + do) cos9^ 

= pdgSffj,, (A.3) 

where 


^277 


^277 


Fi = J 


sin6< Sffj,{9)d9, 


Y2=-J 


cos 9 5J^{9) d9. 


(A.4) 
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We may now define a new operator S such that 


de{S5f)^Rp_de5h. 


(A.5) 


Explicitly, 


{S Sf){ 6 ) = -{F - ysm9)dg6f{6) - + Bq) sin 6 > + Y 2 Pinv {6 + fi'o) cosB . (A. 6 ) 

Clearly, the diagonalization of on the subspace of functions 6 g with vanishing integral is equivalent to 

the diagonalization of S. 

A direct calculation shows that the 3-dimensional subspace V 3 spanned by functions 


fliB) = PinviB + Bo) sinB , f 2(0) = Pinv {0 + Bq) cos B , fsi^) = PinviS + Oq) 

is invariant for S. In particular, the representation of S on such subspace is given by the matrix 


(A.7) 


J 


V 


(^) 

-F 

0 


F 


J 


J\ 1- 


-y 


-y 

0 


(A. 8 ) 




A little algebra shows that the above matrix has one zero eigenvalue which corresponds to the eigenvector 
{—y,0,F), i.e. to the constant function 


Sfid) = Pi„v{B + Bo) {F - y sin B) = 




yZ 


2tt 


(A.9) 


In this case Sg = dgSf = 0. Consequently, the zero eigenvalue does not occur for the operator Rp^^^- The 
other two eigenvalues of the matrix (A.7) are given by 


A± = — ± 




i\ - 


1 - * 


(i+^A-W) 


- F^ + y^ . 


(A.IO) 


of which at least one has positive real part if J > 0. It follows that for J > 0, the stationary solution (A.l I 
is linearly unstable. 


B Algorithm for the Taylor expansion of the quasi-potential close 
to a stationary state 

We discuss briefly in this Appendix how the Taylor expansion of the quasi-potential of the Shinomoto- 
Kuramoto model may be numerically calculated. This will be possible, in principle, at all the orders of the 
Taylor expansion. 

Let us first concentrat e on th e quadratic term In order to explicitly calculate kernel yP{&,■§), we 

and consider the eigenfunctions of operators R)p. 


proceed similarly as in Sec. 
Ho-. 


4.5.3 


and acting in space 


= akUk 


and 


Rn 


„Vk = akVk ■ 


(B.l) 


Assuming ak to be different, we impose the orthogonality relations (4.99). Eigenfunctions Uk form a basis 


of Ho but they are not orthogonal. Similarly for Vk- The kernel of the quadratic term of the quasi-potential 
may be represented as 


^5° (6», d) = ^ Vk {B) vi {d) 

kj. 


(B.2) 
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where ^ki is defined by 


r2'K 


r2Tz 


^ki = 


Uk{9)d9 / ^p'^{9,^)ul{^)d^^ = {uk\‘^ui). 


(B.3) 


The problem is thus reduced to the calculation of Now, fr om Eq. (4.122), we can easily compute the 
matrix elements of in the basis Uk- Indeed, from Eq. (4.122), 


($ ^)kl = (Wfc|$ ui) = 


2kBT 
at + ai 


d9 {deuk){9) p^nv{9) {d0Ui){9). 


(B.4) 


Note, however, that the matrix (4)“^)^; is not the inverse of ^ki because the basis formed by Uk is not 
orthonormal. This problem can be handled with simple linear algebra. We introduce the matrix of scalar 
products 


Pki = {uk\ui) 

whose inverse will be denoted by {P~^)ki- Then, we construct another matrix 


(B.5) 


(B.6) 


with inverse Bki- Finally, 


^kl = y] Pki Bil ■ 

i 


(B.7) 


Once is known, it is straightforward to write the kernel ip^ in the real space using Eq. (4.126). 

We conclude by observing that, once (p^ is known, one could also numerically evaluate the higher order 
kernels (/?” of the Taylor expansion by solving Eq. ( 4.126) usin g the fact that the basis ^Vk is composed of 
eigenstates of operators see Eqs. (4.123) and (B.l). 
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